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We investigate the spectrum of photohadronically produced neutrinos at very high energies (VHE, ^10^* eV) 
in astrophysical sources whose physical properties are constrained by their variability, in particular jets in Active 
Galactic Nuclei (blazars) and Gamma-Ray Bursts (GRBs). We discuss in detail the various competing cooling 
processes for energetic protons, as well as the cooling of pions and muons in the hadronic cascade, which impose 
limits on both the efficiency of neutrino production and the maximum neutrino energy. If the proton acceleration 
process is of the Fermi type, we can derive a model independent upper limit on the neutrino energy from the 
observed properties of any cosmic transient, which depends only on the assumed total energy of the transient. 
For standard energetic constraints, we can rule out major contributions above 10^" eV from current models of 
both blazars and GRBs; and in most models much stronger limits apply in order to produce measurable neutrino 
fluxes. For GRBs, we show that the cooling of pions and muons in the hadronic cascade imposes the strongest 
limit on the neutrino energy, leading to cutoff energies of the electron and muon neutrino spectrum at the source 
differing by about one order of magnitude. We also discuss the relation of maximum cosmic ray energies 
to maximum neutrino energies and fluxes in GRBs, and find that the production of both the highest energy 
cosmic rays and observable neutrino fluxes at the same site can only be realized under extreme conditions; a 
test implication of this joint scenario would be the existence of strong fluxes of GRB correlated muon neutrinos 
up to ultra high energies, >10^^eV. Secondary particle cooling also leads to slightly revised estimates for 
the neutrino fluxes from (non-transient) AGN cores, which are commonly used in estimates for VHE detector 
event rates. Since our approach is quite general we conclude that the detection or non-detection of neutrinos 
above ~10^' eV correlated with blazar flares or GRBs (e.g., with the Pierre Auger Observatory), would provide 
strong evidence against or in favor of current models for cosmic ray acceleration and neutrino production in 
these sources. 
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I. INTRODUCTION 

Neutrino astronomy may provide valuable clues for the understanding of the properties of neutrinos and their interactions 
at energies in the range 10^"*— 10^^ eV, as well as providing qualitatively new information about some of the most interesting 
cosmic objects. This energy range is also of great interest because it can probe the universe at significantly greater distances 
than is possible with known stellar sources (e.g., the Sun, or Supernovae such as 1987a, produce neutrinos in the ~ MeV range 
through nuclear interactions, which would be difficult to detect from more distant sources due to the overwhelming background 
of atmospheric neutrinos from cosmic ray air showers [|l]]). For these reasons, many of the future neutrino telescopes are 
designed for energies > Te V, where the atmospheric background becomes negligible. Among the most promising and ubiquitous 
astrophysical sources of neutrinos at these very high energies are Active Galactic Nuclei and Gamma-Ray Bursts which have 
in common that most of their energetic emission appears in short, distinct flares. The study of the physical processes determining 
the energy spectrum of such neutrino emitting transients is the subject of this paper 

Above TeV energies and up to about 10^^ eV neutrinos are detected predominantly through the Cherenkov effect in large 
volumes of water or ice, using the mass of the earth to capture neutrinos and looking for traces of upward going muons from 
H conversions; above 10^^ eV air scintillation techniques and large air shower arrays become the most efficient to detect 
neutrino induced, deeply penetrating horizontal air showers, where electron neutrinos have the advantage to generate showers 
which are easier to distinguish from the cosmic ray background (see Appendix^ for details and references). The major change 
of detection techniques at about 10^^ eV motivates the distinction between very high energy (VHE) neutrinos at 10^'*— 10^^ eV 
and ultra high energy (UHE) neutrinos at > 10^^ eV. The most obvious source of UHE neutrinos are interactions of ultra-high 
energy cosmic rays with the universal microwave photon background [|^], which predict a diffuse neutrino flux strong enough 
to be detected in air shower experiments More hypothetical is the prediction of UHE neutrinos from processes associated 
with grand unification scale physics, e.g., the annihilation topological defects In the VHE range, the major contribution 
is expected from Active Galactic Nuclei (AGN) [^-[ll|], which are known to emit very high energy gamma-rays. A possible 
contribution of AGN to the UHE neutrino regime is also discussed in connection with the expected event rates of horizontal air 
showers [p^. Also Gamma-Ray Burst (GRB) sources have been proposed as neutrino sources in the MeV— GeV range [|l3|], 
and at very high energies [p4|]. 
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Neutrino production in hadronic models of AGN and GRB is generally attributed to the acceleration of protons in shocks or 
plasma turbulence, known as Fermi acceleration. These energetic protons then interact with soft background photons to produce 
pions (photohadronic pion production). VHE/UHE neutrinos originate in the decay of charged pions, boosted in energy by 
the original proton Lorentz factor, and maybe by an additional Doppler factor due to the bulk motion of the rest frame of the 
relativistic outflows or jets which characterize these sources. The associated decay of neutral pions leads to the production of 
gamma-rays, and it has been claimed that the observed emission up to > 10 TcV in Doppler boosted AGN jets [ [l5| ] is due to 
this process rather than inverse Compton scattering of photons by energetic electrons |]16[]. A relevant issue is that there are 
alternative models of AGN and GRB where the momentum and energy flux of the relativistic jets or outflows are provided by 
or magnetic fields (e.g., [ p^p^ and references therein), rather than protons, and these would be expected to have negligible 
neutrino fluxes. The positive identification of VHE/UHE neutrinos from AGN and GRB would be an indication for baryon 
loaded outflows. The energetic protons may also contribute to the highest energy cosmic ray spectrum, which is observed up 
to 3x10^" cV [|l9|], where AGN and GRBs are considered among the most plausible sources. Here also, the observation or 
non-observation of energetic neutrinos would be a crucial test for these models. 

The details of the photohadronic production of neutrinos via pion decay depend strongly on the properties of the source, 
i.e., its size, lifetime, magnetic field, etc. The magnitude of these quantities are estimated from the observed variability, flux, 
and the measured or inferred distance of these sources. GRB usually last 0.1—100 seconds, but show intrinsic variability 
down to milliseconds, while AGN emit most of their energetic radiation in strong flares lasting several weeks, with intrinsic 
variability on time scales of days down to less than one hour. The transience of energetic emission could improve the association 
of detected neutrinos with their putative sources, because one could use both arrival direction and arrival time information, 
allowing statistically significant statements even for total fluxes below the background level. 

On the other hand, transience and variability sets constraints on the maximum energy of the neutrino spectrum. In the literature 
so far, this has been connected to the maximum proton energy using simple kinematical relations. As we show in this paper, 
however, in these astrophysical scenarios the secondary particles in the photohadronic cascade, i.e., pions and muons, have to 
be considered separately, since cooling processes can have a significant impact on their final distribution. Moreover, one needs 
to evaluate carefully the competing proton energy loss processes that do not lead to neutrinos, which can cause breaks in the 
neutrino spectrum that are not present in the proton spectrum, and thus strongly limit the predicted fluxes at very high energies. 

Starting with a general treatment of photohadronic neutrino production in variable sources, we derive a general upper limit 
for the maximum energy of neutrinos produced by photohadronic interactions of Fermi accelerated protons in cosmic transients, 
which only depends on the total energy of the transient and observational parameters, like duration or (photon) luminosity. 
We then apply our results to hadronic AGN and GRB models, and find that they impose severe constraints on their possible 
contribution to the UHE neutrino spectrum. Notational conventions used throughout the paper and frequently used symbols are 
explained in Appendix^. 



II. PHOTOHADRONIC NEUTRINO PRODUCTION 

A. Proton cooling and neutrino production efficiency 

Photohadronic neutrino production is a result of the decay of charged pions originating from interactions of high energy 
protons with ambient low energy photons. It is accompanied by the production of gamma-rays from neutral pion decay; the 
details of the process are described in Appendix^ for a target photon spectrum following a power law with index a above a 
break energy eb. For the proton gas, pion production acts as a cooling process, and is in competition with other cooling processes 
like Bethe-Heitler pair production, synchrotron radiation, cosmic ray emission, and adiabatic losses due to the expansion of the 
emission region. 

While neutrinos are exclusively produced by charged pion decay, gamma-rays are produced by a variety of processes: besides 
neutral pion decay, tt" — > 77, the major hadronically induced channels are synchrotron radiation from (a) the UHE protons 
themselves, (b) the electrons, muons and charged pions in the photohadronic decay chain (see Appendix]^, and (c) the Bethe- 
Heitler pairs produced in pj pe+e" interactions [pO[]. If synchrotron cooling of secondary particles is negligible, about 25% 
of the energy in charged pions is converted into gamma radiation by synchrotron cooling of the electron produced in muon 
decay [ p"6pl| ]. Normally, the first generation gamma-rays cannot leave the emission region, but rather induce an electromagnetic 
cascade through pair production with low energy background photons, and subsequent synchrotron radiation of electrons and 
positrons. They cascade down in energy, until they eventually escape below some critical energy where the emission region 
becomes optically thin [^]. Hadronically induced gamma-rays are usually in competition with synchrotron and inverse Compton 
photons radiated by primary energetic electrons. 

Cosmic rays can be ejected in essentially two ways: (a) if the emission region has a sharp boundary beyond which the 
magnetic field drops rapidly, protons scattered across the boundary would be ejected, and (b) secondary neutrons produced 
in p7 interactions can escape if (bl) their decay length in the comoving frame is larger than the size of the emission region. 
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CTn ^ R with r„ = 'jpT^^, and (b2) their probability of re-conversion to a proton by, e.g., a reaction nj pTT~ is small, 
expressible by ci„^p ^ i?. Process (a) depends on the detailed structure of the emission region and is usually in competition 
with adiabatic cooling, which affects charged particles due to the adiabatic invariance of the quantity Br? during the Larmor 
motion of the particle = E/eB is the Larmor radius) in a magnetic field decreasing with expansion [g2p. In an isotropically 
expanding emitter w ith con served total magnetic energy this means B (x i?^^, and thus E oc R^^, but other dependences 



may apply (see Sect. fllBJJ). Process (b) is tightly connected to neutrino production, because the dominant channel producing 
charged pions, p7 n7r+, is also the dominant channel for neutron production. The time scale for proton-neutron conversion 
is tp^n ~ tn^p ~ "^tp.TT, which is larger than the time scale for charged pion production, t„± w |^p,7r (Appendix^, because 
there the process p-f ptt+tt" contributes considerably. 

The efficiency of neutrino production depends on (i) which fraction of their energy protons convert into charged pions, and 
(ii) the fraction of energy pions and muons retain until they decay. Condition (ii) can be quantified by introducing the efficiency 
of energy conversion from the originally produced pion into neutrinos, = ^{'^'^° f'^^') and C,^^^ = i(7^'/75J'), f^i" muon and 
pion decay, respectively, where « 7^ is the Lorentz factor of the pion at productio n, and 7^"^ < j!^" < 7J' are the Lorentz 
factors of the muon and the pion, respectively, at their decay (see Appendix^ and Sect. [II C). Similarly we can quantify (i) by 
introducing the charged pion production efficiency, Ctt* = tp/t-K± ^ where tp is the total cooling time of the proton. This leads to 
the total neutrino production efficiency 

1 7^\ 1 tp 



V3 7p 6 7p / ip,^ 

The total proton cooling time is determined by the inverse sum, — - 1^], extending over all participating cooling processes. 
To classify the cooling processes by their dependence on 7p, we introduce total cooling times for photohadronic interactions, 
tp-y, synchrotron radiation, ip,syn, and external cooling processes, tec- Under external cooling processes we subsume adiabatic 
cooling, which has a time scale tad independent of the proton energy, and direct ejection of protons from the emission region. 
The latter may be dependent on the proton energy if diffusive losses are relevant; in the simplest case, however, we can assume 
that protons are confined over the time scale set by adiabatic expansion, i.e., tcsc.p ^ ^ad, which means that 

tec ~ iad = const (7p) . (2) 

The synchrotron loss time can be written as 

tp syn — ^b.svn ( ) with tb,svn — ~. 5 ; (3) 

V7b/ 4rpw2p7b 

where rp — /mp(? m 1.5x10^^^ cm is the classical proton radius, and Ug^p = eB/nipC is the cyclotron frequency of the 
proton. The characteristic proton Lorentz factor used here for normalization, 

7b = ej'if/2eb (4) 

expresses the limit above which all photons in the power law part of the spectrum are boosted above the reaction threshold for 
pion production, assuming that the photon number spectrum can be described as a power law, dNp\i cx e^" de, with an index 
a > 1 above some break energy eb (see Appendix^ for details). The cooling time for pion production can then be written in a 
similar way as 

ip.TT = ib,7r(7p/7b)^"° for 7p < 7b (5a) 

tp,Tr ~ ib,ir(A^ph,b/A^ph) for 7p > 7b , (5b) 

where TVph is the total photon density, A^ph.b is the density of photons with e > Cb- ib,Tr is the pion production cooling time 
for protons with 7p = 7b, and can be expressed as tb.Tr — [cA^ph.bHa]^^, where Hq is the inelasticity weighted effective cross 
section for pion production, as defined in Appendix^ 

The time scales for other photohadronic cooling processes, including neutron ejection, can all be expressed in tp^^. The 



cooling time for the Bethe-Heitler process can be evaluated similarly to Eqs. (g) and (B2); the inelasticity weighted cross section 
for the Bethe-Heitler process is Hq bh ~ Ha/125, for a break Lorentz factor 7b. bh ~ 7b/140, leading to /bh = tp,7r/ip,BH ~ 
exp(5a — 10) for 7p < 7b. bh, and /bh ^ 1 for 7p ^ 7b, bh- To quantify the time scale for energy loss due to free neutron 
escape, icsc,n, we introduce the probabilities for neutron-to-proton reconversion within the length scale of the emission region 
R, due to beta decay, Pnp — exp(— cr„/i?), and due to pion production, Ra-y^p ~ exp(— 2cip.7r/i?), using t„_p ~ 2tp,T^. Since 



the typical energy ratio of the neutron to the pion in a p7 ^ n7r+ reaction is ^4, we finally get icsc.n ~ jPcsc ntp^ir^ where 
-Pcsc.n = (1 - ^'n7— p)(l " Pnp) is the probabifity of the neutron to escape. With t~^ = t^ j^ + 1~ J. „ + t^^^^, and using the fact 
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that beta decay and photohadronic interactions for neutrons are likely in different Lorentz factor regimes, viz., PnpPwy^p 1. 
we can write 

^ 1^ « 3 + exp(5a-10)-2exp(^-^^) -2cxp(^-^) (6) 

for 7p ^ 0.0l7b. For 7^ > 7b, the term expressing the Bethe-Heitler efficiency, exp(5a — 10), is absent. Assuming t^c ~ ^ad 
we can write 

~^ ~ (^P7 ^p,syn ^ec ^^P)"' ~ max(/ad, /p7, /syn) = /max i (7) 

where fsyn = tp.7r/^p,syn and /ad = fp,7r/iad are defined analogously to fp~^, i.e., expressing the energy dissipated in the various 
cooUng channels in units of the energy lost in pion production. The approximation by the maximum-function is best if one 
cooling process clearly dominates, and is useful for the following, qualitative discussion of the spectral distribution of emitted 
neutrinos. 



B. Shape of the time integrated neutrino spectrum 

Because of the low detection efficiency of neutrinos at earth, it is impossible with present techniques to observe short scale 
time variability of cosmic neutrino spectra. Therefore, for an outburst active over a limited time, it is more meaningful to 
calculate the time integrated neutrino count spectrum, rather than the spectral count rate at a fixed time. This also simplifies the 
theoretical treatment, because we do not need to perform a self consistent calculation of the accumulated proton spectrum at a 
specific time — the energy input at any specific energy is simply given by the time integrated proton injection spectrum. 

Clearly, the proton spectrum injected by the acceleration process is not directly observable. We will follow here the scenario 
assumed in most models, that the average injection rate for energetic protons follows a power law in energy, (dNp) = Ipjp'^djp, 
extending from some minimum Lorentz factor 7p to a maximum Lorentz factor jp ^ 7p. We assume that the injection is active 
over a time Ti„j, and that the injection spectrum does not change in time. Then the total, time integrated energy density of 
injected protons is given by 



Up = Tinj / mpC^^p{dNp)^^ = bpiTipC^Tinjlp^l'^ (8) 

where hp is the bolometric correction factor of the proton spectrum relative to its energy content at the highest particle energies, 
given as h,, = \{%,/jpy^^ - l|/|.s - 2| for s^2 and bp = ln(7p/7p) for s = 2. 

Protons injected at a specific Lorentz factor 7p produce charged pions with 7^' « 7p at a rate t~l, over a time tp. Thus, 
the total number density of charged pions produced in the emission process is dN^ = fTlnj (dNp) tpt~l.. Each charged pion 
produces two muon neutrinos {Vf^Vp) and one electron neutrino {v^ or v,,), each with an energy w \'mTrC?^p if we assume 
that 7^° « 7^"= « 7p. Then, the total time integrated neutrino power at the energy Ei, emitted by the source in its rest frame is 

dN^ _ Vup ^/ AE^ 

d\iiE„ nip bp^p~^ \m7rC" 



W)-i^.3^-^r^$h^) , (9) 



where V is the volume of the emission region. The spectral shape is expressed by the function <I'(7p) = 7p *C;/(7p)> which can 
be written as $(7p) = 7p~'*/mil;x(7p) ^ 7p the case of one dominating cooling process; the power law index q is, depending 
on the dominant cooling process, given by 
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(9a) 



where a is the target photon spectral index. For a > 2, an additional spectral modification will occur due to the drop of the 
Bethe-Heitler efficiency between ~ 0.0l7b and ~ 7b, if photohadronic cooling is dominant in this region (/max = fp-y)', in this 
case, one would expect a rapid rise in the neutrino flux in the regime < [30 MeV] 7b. 
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The energy dependence of neutrino event rate as observed in a given detector follows closely the neutrino power spectrum, 
as shown in Appendix^, if we properly account for energy shifts due to Doppler boosting or source redshift, i.e., we expect an 
increase of events with energy for g > 0, a decrease for q < 0, and g w would indicate an event rate almost independent of 
energy. The dominant proton cooling process in the source, and thus the value of q may change with energy. Obviously, such 
a transition of cooling processes at some energy, leading to a spectral break, can only occur if the process taking over has a 
cooling time decreasing faster with energy. Therefore, the spectrum steepens at each break, and for 7p < 7b the only possible 
sequence of break energies is Epj < Esyn(p} ^ Ep for a < 2, and Esyn{p) !i < Ep for a > 2, while only one break can 
exist for a = 2 or above 7b [cf. Eq. (|9a[)]. However, depending on the source properties it may b e that some cooling processes 
are never dominant, so not all possible break energies may appear in the spectrum [cf. Sect. IV A and Fig. |l] for an example]. If 
the proton spectral index is close to the canonical value for shock acceleration, s « 2, we would generally expect an increasing 
event rate at low energies in the regime of adiabatic cooling dominance, and a flat behavior {q « 0) if photohadronic cooling 
becomes dominant, which can be most easily understood as a saturation of the efficiency, Eq. (|l]). The efficiency can decrease 
again if proton synchrotron cooling becomes dominant (for a < 2), unless the proton spectrum cuts off first . The possible 
cooling of secondary particl es in the hadronic cascade may lead to additional breaks, as discussed in Sect. [II C 2 and illustrated 
in an example in Sect. [V B and Fig.0. 



III. NEUTRINO EMISSION FROM COSMIC TRANSIENTS: GENERAL THEORY 



Although neutrino bursts themselves may only be observable in their time integrated appearance, the accompanying burst of 
photons is observable in much greater detail, and allows to constrain the physical parameters determining neutrino production. 
We have to distinguish between low energy photons, which are in general explained by synchrotron radiation of electrons co- 
accelerated with the protons, and high energy gamma radiation, which may be dominated by hadronically induced cascades, 
as discussed above, but could also originate dominantly from inverse Compton photons produced by the electrons. The terms 
"low energy" and "high energy" are here used only in a relative meaning — the absolute energy range for electron synchrotron 
radiation on the one hand, and electron Compton radiation or hadronically induced photons on the other, depends strongly on 
the physical conditions at the source. The following discussion is more focussed on the low energy photon component, which 
is relevant as the target population for photohadronic neutrino production. However, if the high energy photons are of hadronic 
origin, their variability can also give valuable clues on proton cooling times. 

Hereafter we distinguish between physical quantities defined in the comoving frame of the source and observed quantities 
(see Appendix for notational conventions), where we account for a possible boosting of the radiation emitted from the source 
with a Doppler factor T) = [T{1 — f3r cos Ovicw)]^ for a relativistic flow with Lorentz factor F and a velocity /3rC under an 
angle 0view to the direction of the observer. We do not take into account cosmological redshift effects, and just note that they 
might be considered by replacing V = V / (1 + z), if V is the Doppler factor of the emission region in the cosmologically 
comoving frame at the source. The task of this section is to constrain the physical properties of the emission region by observable 
quantities, in order to discuss the various processes limiting the neutrino energy. 



A. Variability time scales and the size of the emission region 

1. The causality limit 

If a flare occurring in a relativistic outflow, boosted with a Doppler factor P, is observed to have a duration T, the time scale 
of the burst in the comoving frame of the fluid is T = TV. This flare time scale covers (a) the time scale for the injection of 
energetic particles, Tjnj, (b) the time scale over which the particles convert their energy in radiation, Tjad, and (c) the crossing 
time the photons need to leave the emission region in the direction of the observer, Tcr- The partial times normally do not simply 
add up, but by order of magnitude the estimate 

T max(Tinj,Trad,T'cr) (10) 

applies. The crossing time is naturally connected to the (comoving) linear size of the emission region; if the emission region is 
not spherically symmetric, we can only limit the comoving size along the line of sight, by the observed duration as 

= cTcr = CprTV , (11) 

where the factor Pt < 1 considers the effect of a delayed emission due to finite injection or radiation time scales; pj. = 1 
means that the emission is homogeneous and instantaneous within the size R — the observed duration is then simply the time 
between the first and the last photon reaching us. Since the condition T > is equivalent to the requirement that the emission 
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throughout the emission region is due to one, causally connected process, Eq. ( |l l| ) may be called the causality limit for the size 
of the transient source. The projected (or lateral) comoving size, R±, is not constrained by the variability time scale, but plays a 
role for the determination of the internal radiation density of the emission region from its observed, isotropized luminosity at a 
specific photon energy, C{e) = Airdle'^ {dNpi^ / de)ohs = 'D'^L{e) with e = Ve, where c^l is the luminosity distance of the source 
(note that the redshift is absorbed in the Doppler factor, as explained above). To account for this, we introduce a geometrical 
eccentricity parameter, by writing the luminosity as 

L{e) = iTTRf^cx^€\dNpi,/de) , (12) 

where e^{dNp\i/de) is the specific energy density of photons with energy e in the rest frame of the flow. For a spherically 
symmetric emission region = 1, while a disklike emission region (iJy ^ i?^) would be describedby a:i'^(_R^/i?ii)^^l. 

The radiation time scale is obviously equivalent to the cooling time scale of the radiating particles; for radiation processes 
involving electrons, it is usually very short. Trad ^ Tci , justifying its neglect. If we consider radiation produced in photohadronic 
interactions, we can write Trad ~ tp, since the time scale over which the electromagnetic cascade evolves can be considered as 
short compared to ip. Tjnj is the time over which an acceleration process is active, e.g., the lifetime of a shock. Obviously, this 
sets a limit on the acceleration time of the particles, t^cc < Jlnj < TV. Also here, this is barely relevant for electrons, but sets 
an important limit for protons. 



2. The internal shock scenario 



As an example how observed time scales in transient emission phenomena may be connected to the size of the emission region, 
we discuss the scenario of energy conversion in relativistic flows by internal shocks. This scenario was suggested originally by 
Rees p3| ] for AGN jets, and was later also applied to Gamma-Ray Bursts [p4[]. 

We consider two plasma blobs of similar mass and density emitted within an unsteady flow at times ti and t2, /S.t ~ t2—ti > Q, 
with respective Lorentz factors Fi and F2, F2/F1 > 1, i.e., the second blob has a larger velocity and thus catches up with the 
first after some time ~ F1F2 Ai. Assuming that their relative velocities are supersonic, two strong shock waves moving in 
opposite directions form when the blobs merge; they are called forward shock and reverse shock, respective to their direction of 
motion relative to the flow. In their center of mass frame (CMF), which has a Lorentz factor F « x/FiT^ in the observers frame, 
the shocked material in region between the two shocks is at rest, and is the source of the radiation. The shocks move each with 
a velocity /3shC, Psh ~ — Fi/I^, corresponding to a internal shock Lorentz factor Fgh ~ \/^2/^i- The linear size of the 
emitter in the direction of the flow, after the merging is complete, is i?,, = 2i?|| / y^, if is the length of the blobs in this direction 
in their respective rest frames, and Xp = 4AF + 3 is the compression factor p^]. Therefore, Ti„j w R'^^/cXpfish = 7cr/2/5sh, 
which is the time each shock needs to cross half this distance. For transrelativistic internal shocks, f3sh ~ the crossing time is 
therefore a good measure for the injection time scale. The total efficiency for the dissipation of energy by the shocks is given by 
^sh = 1 - 2Fsh/(l + F^j^), and is about 20% for F^h « 2. 

We now assume that the radiation time scale is much shorter than the dynamic time scales involved in the shock merging. 
Then, the emission follows closely the motion of the shocks; if the observer is placed at an angle Oview <C 1 to the flow 
direction, the emission of the forward shock appears as a peak of duration 7t = Ti„jP^^(l — (3sh) in the observer frame, while 
the emission of the reverse shock causes a peak with a duration %■ — Ti„jl?^^(l + j3sh) and comparable total energy. For 
Psh ~ the two peaks can thus have different lengths (but of the same order of magnitude), and the total, superposed peak 
might appear asymmetric, with a rise time « Tf/2'D, and a decay time « (7^. — ^T{)/T>. The crossing time is then 
correctly estimated by = TV (or — I), if we define 

T = T-T ^ T,~% . (13) 



We stress that this result is independent of fish- Eq. (13) makes use of two assumptions: (a) that the two plasma blobs have 
comparable densities, and (b) that Trad ^ Tcr. Condition (b) is mostly fulfilled, if we consider synchrotron of inverse radiation 
from energetic electrons. Condition (a) can be assumed to nearly fulfilled in internal shocks — which is the most important 
difference to external shocks, where the densities are usually very different. If (a) is not fulfilled, forward and reverse shock have 
different velocities in the CMF, and also different efficiencies in energy conversion | p^ , p6| ]: hence, the forward and backward 
peak have very different strengths, and the correlation of T and T with % and 7^ is less straightforward. If (b) is not fulfilled, 
i.e., if Trad > Tci, we expect Trad ~ T{V, and Tcr ^ %V. A similar situation arises for fast cooling, but the presence 
of a secondary acceleration process which is not associated to the shock waves (e.g., second order Fermi acceleration, see 
Appendix^, which can keep up a population of energetic particles homogeneously over the region of the shocked gas, and thus 
extend the emission as in the case of slow cooling. Despite these ambiguities, we may assume that for flares with considerable 
asymmetry, T^—% > %^ Eq. ( p^ ) gives a reliable upper limit on the crossing time. We also note that the geometrical eccentricity 
parameter, x^, as introduced in Eq. ([T^, satisfies the relation = X^pX'^, if x',^ is the eccentricity of the blobs before they merge. 
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Since xj, ^ 10 for transrelativistic shocks, this may give rise to assume rather disklike geometries in the internal shock scenario; 
this conclusion, however, may not be over-interpreted because we can obviously not rule out that the blobs have been originally 
elongated in the flow direction, i.e., x'^ < 1. The internal shock mechanism can readily be applied to spherical (or quasi- 
spherical) outbursts, where the up-catching "blobs" have to be replaced by shells emitted with different velocities at diffe rent 
times |24|. We discuss this scenario, and its implication for the geometrical factors Pt and x^, in more detail in Sect. [VB 1. 



B. Maximum energy of accelerated protons 



The predominance of power laws in nonthermal emission spectra suggests that the radiating particles gain their energy by 
a stochastic process. Based on an original idea of Fermi [^, the most commonly discussed stochastic acceleration processes 
fall into two parts: (a) f irst order Fermi-acceleration by diffusive scattering of particles across strong shock waves, also called 
shock-acceleration p^ ] and (b) second order Fermi-acceleration, where the particles gain energy from the scattering at plasma 
waves 1 29 , 3^. Since plasma waves are responsible for the scattering, and thus for the diffusive motion of of particles in shock 
acceleration as well, it is most likely that both processes combine if strong shock waves are present [^]. Fermi acceleration 
is assumed to be the dominant energy dissipation mechanism in AGN cores and jets, and in Gamma-Ray Bursts. Since Fermi 
acceleration works independent of particle mass and charge, any protons or ions present at the shock should be accelerated as well 
as electrons. It has been claimed for various classes of objects that this could be the generating process of the observed cosmic ray 
spectrum, up to the highest energies of order 10^*^ eV [[l9|]. Here we discuss the maximum energies of protons Fermi-accelerated 
in emission regions constrained by variability. The time scale for Fermi acceleration is expressed as a multiple of the Larmor 
time of the particle, iacc = ^f^l — 27r6'F?'L/c, where we assume 9p as constant for simplicity; as discussed in Appendix^ this is 
only true for special assumptions on the diffusion coefficient (e.g., Bohm diffusion). Concerning its magnitude, 6^ ^ 1 applies 
in the most cases, but 6^ ^ lis probably possible for acceleration at relativistic shocks (Appendix^. 



1. Larmor radius and adiabatic limits 

To be accelerated up to an energy Ep by a Fermi mechanism, we have to require that the protons can be magnetically confined 
in the emission region, viz., Ep < ei3i?niin, ^min = min(_R|| , i?^). Using the relations i?,, < cTV, Ey < jiriTrc'^jp, and 
Si' = E^T), we obtain the limit 

£. < 4. = ^ BTV^p^p, . (14) 

The factor Pl < 1 considers that usually only a limited fraction of the effective size i?niin of the emission region can practically 
be used for particle gyration; we will assume ^ ^ in the following. 

Similar limits can be derived from acceleration time constraints arising from the dynamic time scales involved in the ac- 
celeration process. The first condition of this kind is iacc < which can be specified in the internal shock scenario as 

< Tci/29-p(3sh, and is obviously equivalent to Eq. ( p4| ) with /5l — (26'F/3sh)^- The second condition is the limitation of the 
acceleration time by adiabatic cooling in a decreasing magnetic field, t^cc < ^ad = 2|i3/-B|. In an expanding emission region, 
we usually find B cx R^", where a > and R is some characteristic size of the emission region. In the general case, in 
particular for non-isotropic expansion, a may depend on the choice of R; in an isotropically expanding emitter with conserved 
magnetic energy we have a = 2, which we may use as a canonical assumption hereafter Defining the velocity of expansion 
as /3cx = R/c, this results again in Eq. (0), with Pl ~ {'^ct9p/3cx)^^ ■ Adiabatic cooling is most relevant in a freely expanding 
relativistic fluid, /3cx ~ 1, or for rapidly decaying magnetic fields, a ^ 1 (however, we have to assume that the field decay is 
adiabatic, i.e., \B/B\ ^ For second order Fermi acceleration, both constraints, tacc < 21nj and iacc < ^ad, are equivalent, 
because the injection time is limited by the adiabatic drop of the Alven speed, which leads to Ti^j ~ tad- 

In conclusion, Eq. ( p4| ) with Pl ^ min(i, 0^^) applies for our canonical assumption that the involved hydrodynamic pro- 
cess are at least transrelativistic. Acceleration time constraints then dominate for 9p ^ 1, which is particularly the case for 
second order Fermi acceleration, or acceleration at nonrelativistic, quasi-parallel shocks (cf. Appendix^, while for any faster 
acceleration mechanism the geometrical extension of the emission region, constrained by variability, sets the limit on 



2. Limits due to radiative cooling 

Synchrotron cooling of the protons during acceleration limits their maximum energy through the condition tacc < 
Writing the acceleration time as tacc = 2tt9-p^p/ul)b,p and using Eq. (^, we find 7p < 3/ y^Svr^^J^ryl^^, leading to 
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3 m^c^D I e 

< f.,sy„(p) - 8^;~V2^ ■ ^^^^ 

In the same way, proton acceleration must be faster than the coohng due to photohadronic interactions, iacc < ~ fp^tp^-^. 
To express tb,7r by observable quantities, we use Eq. ( p^ to write 

ib,7r = 7^ With Tc = — —2 , (16) 

where £b = '^(eb) is the isotropic luminosity of the source at the observed spectral break, related to the break energy in the 
comoving frame by eb = ^h'D^ and H2 ~ 22 /ibarn (cf. Appendix^. Using Eqs. (^ we obtain from tacc < tp-y 



Of Spy -Lc 



^'pTcT'b 



where we inserted the Doppler scaled characteristic Lorentz factor, 7b = ef^/2eh, which does not depend on V and is related to 
the comoving characteristic Lorentz factor by 7b = 'jhT^, and 

r - j T^fpj/2TpT for % < 7b 

^ 1 T^/p7Afph/2iVph,bTp^ for % » 7b ' ^ ^ 

The case of Eq. ( p^ obviously corresponds to % < 7b ; the case % » 7b is described by setting a = 1 and using the proper value 
for [note: in Eq. (|l6|), the actual power law index a has to be used in any case]. It should be noted that Eq. ( [l7| ) only considers 
the photon density in the burst connected to its intrinsic luminosity. If, however, the relativistically moving "blob" is embedded 
in an ambient photon field, which is isotropic with respect to the observer, the photon density seen in the comoving frame of 
the blob can be considerably higher than inferred from the observed luminosity by Eq. (|l2|). The reason is, that this additional 
component would appear unboosted for the observer, and therefore probably only as a small fraction of the apparent luminosity 
of the emitting blob which is boosted by a factor V^, while in the comoving frame the photon number density of the ambient 
component is increased by a factor T) due to Lorentz contraction and may thus dominate the photon density. This scenario might 
be relevant in AGN, if the emission region in the jet is close to the AGN core, and if the core radiation is isotropized by a plasma 
halo (see, e.g., [ B2p), and photohadronic interactions may limit the neutrino energy to values significantly below the upper limit 
expressed by Eq. (|17[). 



C. Cooling of secondary particles in the hadronic cascade 

1. Cooling processes and time scales 

Pions and muons are weakly decaying particles, with comparatively long lifetimes, t^^ = 2.6x10" , ^ 

respectively. For secondary particle Lorentz factors 7^' ~ 7p 10^~^, which are quite reasonable and readily considered in 
most models, their lifetime in the comoving frame of the emission region, = ''^^''7*, can be of the order the dynamical time 
scale of the flare (e.g., in Gamma-Ray Bursts). Moreover, their synchrotron losses are by a factor [mp/miJf' ^ 10^ stronger than 
for protons. Adiabatic cooling of muons has been considered for neutrino emission of Gamma-Ray Bursts [|l4||, and synchrotron 
cooling of pions and muons has been discussed for extremely magnetized environments, e.g., neutron star magnetospheres [|33|]. 
Most of the literature about neutrino emission of AGN, however, neglects this effect. We will show here that this is not justified, 
and derive the critical Lorentz factors, 7^, above which the energy loss of muons and pions plays a role. Obviously, we have to 
distinguish between neutrinos from the decay of pions and muons, because of their very different lifetimes. Energy losses of the 
muons are generally more relevant, which affects in particular the electron neutrinos, arising exclusively from muon decay. This 
is important for neutrino detection in UHE air shower experiments, where electron neutrino showers are easier to distinguish 
from the atmospheric background and are therefore proposed as providing most of the expected signal [ p"2p4| ]. 

Secondary particles cool adiabatically prior to their decay if r,t > 2\B/I3\, which gives a critical Lorentz factor 7^, — 
2i?/a/3oxrr for S oc leading to 

where is the neutrino energy in the rest frame of the decaying particle, which is roughly |;m^c^ w 35MeV for pion 
decay, but slightly smaller (BOMeV) for muon decay and may therefore be distinguished (see Appendix Bh. Analogously, 
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pions and muons can undergo efficient synchrotron coohng if < t^.syn- The critical Lorentz factor is thus found from 

7b^b,syn(™*/TOp)^ = 7*'^*''' yielding 



S^.syni*) ^ — EZ-D with = /__ , (20) 



where = e^/m^c^ is the classical radius of the particle, and Wi,.* = eB/rrii^c its cyclotron frequency. The characteristic 
frequency, w^, for synchrotron losses of pions and muon is found as oj^r = 5. Ox 10^^ s^^ and u)^ — 4.7x 10^^ s^^, respectively. 
Secondary particles may also suffer inverse Compton (IC) losses from interactions with background photons. The corresponding 
cooling time is related to the synchrotron cooling time by the well known relation 

t - t ^ - t ^'^'^^P'' (21) 

where 6ph is the bolometric correction factor relating the total photon luminosity over the spectral range Cb to e to Lb given by 
&ph = K^/^b)^ " — 1|/|2 — a\ for a 7^ 2, and 6ph ~ ln(e/eb) for a = 2. For high pion and muon energies Eq. (21 ) might be 
modified by Klein-Nishina corrections, leading to a suppression of IC cooling. 

Unlike the proton case, IC cooling is the most relevant process for photo-interactions of pions and muons. Their IC cooling 
time is related to the proton IC cooling time as ic = (rrhi^/mpYtp^ic, while the Bethe-Heitler cooling time only scales with 
i,t,BH = {'mi,/mp)tp^BH- For a — 2, one can show that i*^BH ^ St^^ic, if Klein-Nishina corrections are disregarded. For the pion, 
there are additional channels due to meson resonance excitation; the lowest energy process is tt^j ^ ^ tt^tt", which has a 
theoreticalpeak cross section of ~ 50 /zbarn (determined from the — > 7r*7 decay branching ratio by use of the Breit-Wigner 
formula 03]), at a photon energy of e^^p ~ 2 GcV in the pion rest frame. Compared to the pion production off the proton 
via the Z\-resonance, where the characteristic photon energy is ep^A ~ 300 McV, this interaction is suppressed by a factor 
^ ^(£p^zi/e7r^p)°; for a w 2, and neglecting the finite hfetime of the pion, secondary photon scattering has a probability of 
< 0.3%, and can thus be neglected. However, the process may be relevant in inverted photon spectra, e.g., for pion production 
in the Rayleigh-Jeans part of the thermal background. 



2. Spectral modification at the critical energies 



The "maximum" neutrino energies derived in Sect. Ill B arise from the balance of energy gain and loss processes. The 



stochastic nature of these processes allows the particles to exceed such "limits" with some, usually exponentially decreasing 



probability. This has been shown for Fermi accelerated particles subject to synchrotron losses |36|, which show a largely 
unmodified extension of the power law spectrum up to the cutoff energy (defined by balance of gains and losses), followed by an 
exponential-like cutoff; depending on the detailed parameters, a pile-up may occur at the cutoff energy. Although the stochastic 
behavior of photohadronic losses is quite different from synchrotron losses, we may expect a similar result in this case, and 
also for adiabatic losses. In general, we can assume that if the neutrino energy is limited by the maximum proton energy, the 
spectrum will continue as an approximate power law up to £,y, and drop off rapidly for E^, > 

The situation is different for the decay of unstable particles, where the particle can decay within a time At after its production 
with a probability Pdc = 1 — exp(— Ai/r*). Far above the critical Lorentz factor, 7^° ^ the decay probability within a 
cooling time scale in the fluid frame. At ~ icooi(7^) ^ t*, can be approximated as Pdc(^cooi) ~ ^cooi (7* )/''*■ The critical 
Lorentz factor is defined by the condition that tcooi(7*) — Ji^TZ^ which allows us to write tcooi(7*) = ''*(7*/7*)^'"- For 
adiabatic cooling, <cooi = const(7*) then means w — 1, while for synchrotron cooling we have w — 2; cooling by secondary 
photon scattering — if relevant — would correspond to w = a. For pions or muons produced with 7^' < 7*, neutrinos are 



produced with an energy S^, « EZlp'D in the observers frame, and their power spectrum is Ci, oc S^, as discussed in Sect. [IB 
For 7P' ^ 7^, the power spectrum is modified by the probability to decay within a cooling time icooi(7r)' ^i^-' 

(X Pdc(icooi) (X for » EZ%V . (22) 

Therefore, the critical Lorentz factor marks a spectral break of magnitude Aq — —w, rather than an exponential cutoff. Clearly, 
this simplified analytical estimate does not treat exactly the energy evolution of the pions and muons, thus neglects particle 
number conservation. Considering this would lead to a pile-up of decaying pions and muons around their respective critical 
Lorentz factor, before the spectrum turns over into a f^^™ behavior. The strength of the pile-up is correlated to the magnitude 
of the spectral break, and is therefore expected to be stronger for synchrotron cooling breaks than for adiabatic cooling breaks. 
The break due to adiabatic cooling of secondary particles is comparable to the spectral breaks occurring at the transition between 
different dominant proton cooling processes, Eq. (Ba). The Aq = —2 break caused by pion or muon synchrotron cooling. 
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however, is stronger and can be easily confused with an exponential steepening: it causes a drop of events of one order of 
magnitude over half a decade in energy, similar to the exponential function around its critical energy. In practice, we may 
therefore consider £i/,syn(*) as a cutof f energy of the neutrino spectrum and compare it with with the cutoff energies due to 
proton cooling as derived in Sect. Ill B| . 



D. Model-independent discussion of spectral shapes, maximum energies and fluxes 



1. The parameter space 



The free parameters describing a transient fall into two classes: (i) observable parameters, i.e., the characteristic time scale, 
T, and the isotropized luminosity, £, and (ii) theoretical parameters, like B, T), 9p, et c. We consider the former as given 
for any specific transient (disregarding possible problems in their determination, cf. Sect. |^11 A2 ), while the latter can only be 
constrained through general physical considerations or additional observations within a certain range. Eqs. (p4|), (|l5|), and ( |l7| ) 
show that the maximum energies depend strongly on some of these parameters, in particular on the magnetic field B and the 
Doppler factor V. In contrast, the hydrodynamic parameters of the flow, i.e., a, /3cx and /Jgh are of order unity, and can generally 
be recast into some reasonable assumption for Pl- A special role is played by 9p, which describes the speed of the acceleration 
process: as shown in Appendix^ first and second order Fermi acceleration is limited to 9^ > 1, but we could easily consider 
any faster acceleration process in our analysis by inserting the appropriate 9^. We therefore shall use 9p as a fixed parameter, 
with the canonical assumption that ~ 1, while B and T) span the two-dimensional parameter space describing a transient. 
From Eqs. and we can immediately derive a qualitative division of the parameter space: 

(1) For any given V, synchrotron losses of protons dominate over both photohadronic or adiabatic (Larmor) limits for mag- 
netic fields larger than some value i3sy„(p) (P). Analogously, there is a magnetic field i3syn(^)(2?) above which muon or pion 
synchrotron losses dominate over all other proton loss processes; the relation of i3sy,j(p)(I?) and i3syn(*)(2^) depends on the 
other parameters. 

(2) For any given B, photohadronic interactions of protons dominate over both synchrotron and adiabatic (Larmor) limits for 
Doppler factors smaller than some value Vp^ (B). Similarly, there is a limit 'Dpy^^, below which photohadronic interactions also 
dominate over pion or muon synchrotron cooling. 

(3) For any given V, there is a limiting magnetic field i3ad(p) (2?) below which adiabatic or Larmor limits dominate over 
synchrotron or photohadronic cooling of protons, and another value i3ad(*) below which it dominates over either muon or pion 
synchrotron cooling. 

It is much more difficult to determine under which conditions adiabatic cooling of secondary particles dominates; we will see 
in Sect. IV B that, if at all, this can happen only in a very limited region of the parameter space. For a more illustrative discussion 



of the parameter space constraints specific to AGN jets and Gamma-Ray Bursts see Sect. IV A 2 and Sect. [V B 2, and the figures 
shown there. 

To get a quantitative idea about the above parameter space division, we derive the condition under which the delimiting 
energies £^,syn(p), ^v.pi, and E^^i^ are all equal. This is equivalent to the condition iacc = tp.syn = tp^ — 2'Kpi^9pR/c, 
representing three equations which can be solved for the three variables B, T) and 7p; note that the latter term reduces to 
2R/ca(3cx = iad in the case Pl = (Tra/^ox^'F)"^, i-e., that adiabatic cooling rather than space limitations determines the Larmor 
limit. The corresponding solutions for the magnetic field, Doppler factor, and maximum proton Lorentz factor are called B*, 
V* , and 7*, respectively, and are for the case 7p < 7b given as 



V* 

B* 
7p 



3'Y^a — 1 a-\-2 nA — a 



= By 



a+3 ^a+4/i3 



-2y 



7b T^T^+-V^+40-("+')l'' 



(23a) 
(23b) 
(23c) 



with y — (2a + 10) ^, where we introduced the magnetic field scale Bh 
dimensionless quantity 



dmlc-^/Sne^jl = [7.3x lO^i G]7^7^ and the 



9 cT Pt 
Svr Tp^l 



(24) 



The relations for the case 7p ^ 7b are obtained from Eqs. ( p3[ ) by setting a = 1 and use the proper value of from Eq. (p^. 
The corresponding maximum neutrino energy is obviously given as £* ~ ir7i^c^I?*7*. Eqs. ( ^ ) define a unique reference point 
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in the parameter space of magnetic field and Doppler factor, which allows the discussion of the relation of cooling processes of 
the proton independent of any other physical properties of the transient; for notational simplicity, we refer to it as the star-point 
hereafter 

For some applications it is useful to normalize the magnetic field energy density by the comoving photon energy density. We 
introduce the equipartition parameter 

?S7 = - 7^7 , (^J) 

where 6ph is the bolometric correction factor of the photon spectrum as defined in Eq. (|T]). Without protons, ^b^j ^ 1 would 
correspond to approximate energy equipartition of magnetic field and relativistic particles, Ug ^ Uf. ^ Uph, because of the high 
radiative efficiency of electrons. In baryon loaded flows this is not generally the case because of the contribution of relativistic 
protons with low radiative efficiency, and Up + Ue ^ Ub may imply ^b-^ 3> 1 if the acceleration process works more efficiently 
for protons. The corresponding star-point value is 



7'Y^3a+l a — 2/)6— 3a 



2y 



(26) 



2. General upper limits on the neutrino energy 



For any given Doppler factor, the highest neutrino energy can be achieved for B ~ i?syn(p)(I?), because with increasing B, 
^i/^L increases and £i,^p^ remains constant, while £y,syn 
T> <T>* and £,y^p-y = £v,syn(p) for T> >!)*, leading to 



(p) decreases. The equations determining i3syn(p) are £^^i^ = £y^sya{p) for 



syn(p) 
syn(p) 



B*(V/V*) 



-fe 



for V < V* 



= B*{V/V*)-^'^ for V > V* 



(27a) 
(27b) 



with k = 2 + 4/(a + 2) in the case 7p < 7b, and fc = ^ for 7p ^ 7b. If we use ^s-y as a coordinate of the parameter space 
instead of B, and the maximum neutrino energy is reached for CB7,syn(p) — for V > V*, and CB7,syn(p) — 

£^*B^{T) /'D* Y'^'^^ otherwise. This provides an upper limit for the maximum neutrino energy as a function of V 



£v < LAT^) = £t{'D/V*)^+''/^ for V <V* 
£v < £u,p{T^) = £t{VlV*f/^ for V >V* 



(28a) 
(28b) 



Stronger constraints on 5^ may exist from secondary particle cooling. We will only consider synchrotron cooling, bec ause the 
impact of adiabatic cooling of secondary particles on the neutrino spectral shape is relatively weak, as argued in Sect. Ill C 2 ; 
moreover, it is often dominated by the other cooling processes, as we will see in Sect. IV 
for T) < T)*, and £u^sya{*) = ^v,i. otherwise, we find the solutions for the magnetic field 



From the conditions £^ 



syn(*) 



-1/2 



for 
for 



V < V* 

V > V* 



(29a) 
(29b) 



where fc^ = 1 + 3/(a + 1), and B*^^^^^^ = ^ B-rBi,/ piJ)* with Br = mpc/eprT and = m^^w^c/e. Using as a variable, 
we obtain ^b, ~ es7.sy„(.)(2'/2'* ^nd ^b, - ?B^,sy„(.)(2?/2?*)', respectively, with 



S7,syn(*) 



C^plT^Xr.BrB^T>*^ 
2pL^ph^b 



(30) 



The transition from photohadronic to adiabatic dominance at _B = -Bsvnf* )(P) is exactly atV = V* in the case 7p ^ 7b. For 



7p < 7b the exact transition value could be easily found by setting Eq. (29a) and Eq. (29b) equal and solving for T), but the 
approximate division at V* will usually be adequate in practice. Secondary particle cooling determines the maximum energy in 
the case -Bsyn(*) (2^) < ^syn(p) (2^)' and the corresponding upper limits to the neutrino energy are 
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L < LAT^) = £uAT^*){VlV*f'^ for V>V* 



with 



(31a) 
(31b) 



(31c) 



Obviously, V* can be replaced by any reference Doppler factor in Eq. (31b). A more illustrative discussion of the relation 
between i?syn(*) {P) and i?syn(p) (^) in the different regions of the parameter space is given in Sect. [V 



3. Efficiency considerations and neutrino flux limits 

Relevant for the observability of ultra-high energy neutrinos is not only the maximum neutrino energy a specific class of 
sources can provide, but also which total neutrino luminosity is associated to this energy. This requires a discussion of the 
neutrino efficiency as a function of the parameter space. 

If we assume that B < i3sy„(,t), i.e., we are in the part of the parameter space where secondary particle cooling plays no 
role, then the neutrino efficiency can be written as = ^tp/tp.Tr, which can be rewritten as = ^/"^x if '^^^ cooling process 
clearly dominates. In the star-point, where all cooling times are equal at the maximum proton energy, the neutrino efficiency at 
this energy is most easily derived as Q, = (6/^^)^^, with f*^ w 3 + exp(5a — 10) since both terms for neutron reconversion 
in Eq. (^ are small in the star-point. Using the scaling properties of the cooling times, we can then derive the following simple 
expressions: 

If photohadronic cooling dominates, 

Cu ~ ^/p^^min with fpj^min « 1 + exp(5a - 10) , (32a) 

which leads to C^.ad ~ 5 if a < 2, i.e., if Bethe-Heitler losses can be neglected. In general, the total hadronic radiative efficiency 
in this region is w 1, which approximately equal parts (for a < 2) radiated in neutrino and electromagnetic channels. Cosmic ray 
ejection is generally suppressed by neutron reabsorption, and direct proton ejection can usually assumed to be marginal, except 
for special geometries and, maybe, at the highest energies. 
If adiabatic cooling dominates. 



V 



-(3+a) 



c^-'Ci^^j llJ ' ''''' 

where Q is the efficiency at the star-point. Here, most of the energy is not radiated, but reconverted into kinetic energy of 
expansion. Thus, the total hadronic radiative energy decreases in total, while the distribution of the radiated energy between 
cosmic rays, neutrinos and photons remains constant, and approximately 2:1:1 for a < 2. 
If synchrotron cooling dominates, 

— (i)(fr(^r(ir ■ 

Like in the case of photohadronic dominance, the total hadronic radiative efficiency is close to 1, but the predominant part 
of the energy is emitted in electromagnetic radiation. The correction factor 1 < fp-yl fpi < 3 considers that in the region of 
synchrotron dominance neutron reabsorption may or may not play a role, depending on the ratio of the photohadronic interaction 
time scale to the crossing time. 

Note that, instead of the star-point, any point in or at the border of the respective dominance region could be used as a 
reference point in Eqs. (^2|). If secondary particle cooling plays a role, the situation is similar: For synchrotron cooling, the 
total hadronic radiative efficiency remains constant, but more energy is channeled into energetic photons; for adiabatic cooling, 
additional energy is reconverted in bulk kinetic energy. We do not discuss the scaling properties for these cases; rather, we will 
use £u,syn{*) as limiting energy and discuss the efficiencies only for £y < £,y,syn(*)^ and neglect the effect of adiabatic cooling 
of secondary particles altogether. 

The neutrino efficiency alone does not allow to derive flux rates; we also have to make assumptions about the energy density 
of relativistic protons in the source. Here we can use the standard equipartition argument discussed above, and introduce a 
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parameter = Up/ug ~ 1, which allows us to express the time integrated proton injection energy density, Up by the parameter 
space variable ^b-i defined in Eq. (^5|). The neutrino luminosity at the maximum energy can then be expressed relative to the 
bolometric photon luminosity as = Cv /hi,\iChT , and we define 



£„ = £„(X>) 



(33) 



7p = ^pCD) 



where 7^ denotes the maximum proton Lorentz factor attainable for a given T) for which secondary particle cooling can yet be 
neglected. If ^syn(,t) {T^) > -^syn(p) (^)' £.B-y and Ci/ {%) have to be evaluated at i? = -Bsyn(p) C^) and we obtain 



f — 

/ — 

while for V = D* we obviously have £^,p 
be used to determine 7^(1?), leading to 



4&p/p-y, mill \ 



V 



2a-4 
a+2 



4-2a 



^bpfj 



B7 j ^ 
V* 



PI 



for V < V* 



for V > V* , 



(34a) 
(34b) 



^pB£.*By/Qbpf;^. If 4yn(*)(25) < 4yn(p)(2'), ^Bj = ^B-y ,syn{*){'D) haS tO 



^pB^Bj,syn{*) 


(- 






^pB^Bj,syn{*) 


(- 







a+1 



3 — g 
2 



for V < V* 



for V > V* 



(35a) 
(35b) 



We note that, if secondary particle cooling plays no role, £^ is independent of T) for a — 2, and has a maximum for T) ~ V* 
for a > 2. In case of the dominance of secondary particle cooling, a maximum is only obtained for a > 3, while otherwise 
£y continues to rise with V. Therefore, for given C^, both the maximum neutrino energy and the power flux at this energy 
increases with V in most cases, and a general upper limit cannot be stated. However, we have to consider that increasing V, 
while maintaining jp = 7p, implies that also S^s-y and therefore the total energy dissipated into relativistic particles and magnetic 
fields by the transient, (£t, increases. Since €t is usually given, or at least limited by fundamental principles for any specific 
kind of source model, we can impose an upper limit 



< CS7,: 



1 



4n €t 



1 



(36) 



where is the total angle over which the energy is emitted; the case of emission in a thin, freely expanding jet corresponds to 
n w V^^. This can be transformed into an upper limit on the Doppler factor for which both = and C^iSiy) = iubphChT 
can be attained, which we call V; for simplicity, we confine the discussion to the case ^57. max > ^B-y syn(*)' which means 
V > T>*, and assume ^57, max ^ 1- From Eq. (31b) we then obtain a general, model-independent upper limit on the neutrino 
energy as 



n*bp^CbT{l+^pB)eBy 
47r€T 



n*bp^Ci,T{l+^pB)eBy 



syn(*) 



for 4yn(.)(I>) > 4yn(p)(^) (37a) 
for 4y„w(I>) < 4y„(p)(I>) , (37b) 



I and — for ft* = Q = const > V 



where f is givenby Eq. (31c). The power law indices in Eqs. (^7[) are: Kp 

while Kp = Ki, ^ ^ with ft* = V*^^ in the case of a free jet, Q. ~ V^^ . Eq. ( ^7| ) is indeed a true upper limit for 5^: increasing 
T) beyond T) while keeping S^By = £,B-y,niax, implies B cx T>~^ for constant 17, and B cx for a free jet; since the maximum 
energy is determined by the adiabatic limit, Eq. (p4|), this leads to £i, cx and E^, cx const, respectively. It is obvious that any 
other choice of parameters, e.g., P < 2? or £^B-y < $S7.max, leads to lower limiting energies. From Eq. (35b) we obtain 



13 



B-y 



47r£T 



B7,syn(*) 



'^bpfpj 



for i?sy„(*)(X') > Bsy„(p)(X>) 



(38a) 



for 4y„w(I>) < 4y„(p)(I>) , (38b) 



with Ap = (2 - a)/7 and A* = (3 - a)/10 for = r2* = const, and Ap = (2 - a)/4, A^ = (3 - a)/6, and fi* = P*"^ in 



the case of a free jet. Eq. ( P3[ ) is not an upper limit, because higher values are generally allowed for lower Doppler factors; thus, 
one can increase the emitted neutrino power compared to the value in Eq. (|8|) on the expense of the maximum neutrino energy. 
We also emphasize the role of the baryonic energy content: increasing S^pB increases {^{V), but decreases V and therefore 
the upper Umit on We note that the dependence on the value of is low, so that in many cases £* and £* will give good 
order-of-magnitude estimates on the possible ultra-high energy neutrino of a transient. We illustrate this discussion in some 
more specific apphcations in the next section. 



IV. NEUTRINO EMISSION FROM UNSTEADY ASTROPHYSICAL SOURCES 

In this section, we apply the general theory developed in Sect. |l| to specific models of astrophysical transients, which are 
commonly discussed in the literature. Our aim is to explore the parameter space of these models more extensively than usually 
done in the literature, and to check the results for consistency with the limits set by secondary particle cooling, which is disre- 
garded in most papers. To simplify the discussion, we ignore in the following the geometric parameters Xj^ and p-r introduced in 
the last section, and assume = pr = I- However, since all times are normalized by the size of the emission region, and all 
luminosities by the photon energy density in the comoving frame, they can easily be re-introduced by replacing T T p^ and 
C ^ C/x^'in all equations of this section. 



A. Neutrinos from AGN jets 



There are two classes of AGN models which predict neutrino emission, both involving normal e-p flows. One class assumes 
particle acceleration at shocks in the accretion flow very close to the black hole, and produce neutrinos by both p7 and pp 
interactions ||^,p[ p^p7| ] . The other applies to radio loud AGN, which show extended radio jets, and locates the emission region 
at internal shocks in the relativistic jets at larger distances from from the black hole The highest energy neutrinos would 

then be expected from blazars, which are AGN jets pointing in the direction of the observer, because the energy is boosted by the 
Doppler factor V ^ Fjct ^ 1 [|8|. We discuss this class of models, i.e., the AGA^ jet models, in the following; some interesting 
implications of our results on the other class, the AGN core models, are described in Appendix eI 



1. The "proton blazar" scenario 

Blazars are known to emit electromagnetic radiation from radio wavelengths up to the TeV-gamma-ray regime. Their spectrum 
shows a typical "two-hump" structure, where characteristic photon energies depend on the source luminosity: In high luminosity 
blazars, such as 3C 279, the lower hump cuts off at optical wavelengths, while the high energy emission extends up to at least 
10 GeV; in low luminosity sources, like in the nearby objects Mkn421 and MknSOl, the lower hump extends in flares up to 
10—100 keV |39,4^, and the high energy emission is observed up to > 10 TeV 1 15 4l]]. While there is agreement that the low 



energy hump is due to synchrotron emission of energetic electrons, the origin of the high energy emission is unclear: It can 
be explained (a) by inverse-Compton emission of the same electron population producing the low energy synchrotron radiation 
(e.g., [ jl7| , p^ , and references therein), which could arise also if the jets are leptonic (consisting of and magnetic fields, 
with few or no protons); or (b) by electromagnetic cascades induced from the decay of photohadronic pions [|T]]. The latter 
mechanism is referred to as the "proton blazar" or simply "hadronic" scenario, and gives rise to considerable neutrino fluxes 
[^,|l^, while the leptonic models obviously do not. On the basis of gamma-ray observations alone, the issue of the dominant 
radiation process in AGN jets is not settled yet [|43|]; the observation of correlated neutrinos could resolve this issue. 

The typical bulk Lorentz factors Fjet of AGN jets can be estimated from the apparent superluminal motion of blobs in the jet; 
a recent investigation of 43 AGN indicates that Fjct ^ 30 [Q], and the typical inclination angle of the blazar jets to the line of 
sight is inferred to be (6viow)bi ^ 5°, confirming the estimate V ^ Fjot ^ 10 obtained from AGN unification models [Bsh; we 
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will use I? < 30 as an upper limit estimate. A lower limit on the Doppler factor of TeV blazars can be found from the observed 
emission of photons with e ^ 1 TeV, where the emission region must be optically thin with respect to 77 ^ e+e" reactions of 
gamma-rays on intrinsic soft photons, which leads toV > IH^JJ^^^^^^ , or roughly I? > 3. A similar limit is obtained for the 
high luminosity blazar 3C279, using EGRET observations and assuming the emission to be optically thin at 1 GeV [^. 

In contrast to the bulk Lorentz factor, the magnetic field strength in blazars is more difficult to estimate, although hadronic 
models typically invoke i? > lOG based on equipartition arguments [Q. A test of this claim is possible by observing the 
synchrotron-self absorption frequency of the variable emission of blazars; an analysis of the spectral shape and the multifre- 
quency variability of blazars suggests that this turnover is at an observed frequency ~300 GHz |47|. Boosting into the comoving 
frame of the blob with I? < 10, this would be consistent with synchrotron-self absorption for i? lOG and a relativistic 
electron energy density ttg ^ B^/Stt, which is expected in hadronic models ijl^]. 

The emission from blazars is strongly variable, with activity periods taking turns with quiescent periods on a typical time 
scales of months [^sjl; this kind of long-term variability is observed at all frequencies, and appears to be largely correlated. Also 
within an activity period, the TeV emission from Mkn421 and Mkn501 shows clearly separated flares with doubling time scales 
from days down to less than one hour pl| , p9| ], viz., T ^ 10'^— 10^ s, with correlated variability of the synchrotron emission at 
X-ray energies [^0[]. This suggests an identification of these short-term flares with transient, causally disconnected acceleration 
regions of energetic particles, e.g., as expected in the scenario of internal shocks in the jets. 



2. The parameter space for time-integrated neutrino spectra from blazar flares 

We assume that the relevant target photons for photohadronic pion production are the synchrotron photons in the low energy 
spectral hump produced by accelerated electrons, since the number density of photons in the high energy hump is too low and 
can be neglected. We confine the discussion to low luminosity TeV blazars, where the low energy photon spectrum extends to 
e > 1 keV; for high luminosity blazars, like 3C 279, in which the synchrotron component cuts off at optical frequencies, our 
power law approximation is not applicable and a more detailed calculation would be required. The target photon spectrum in the 
comoving frame of the relativistic flow can then be approximated by a power law with a typical index (a) w 1 .7, if we use a power 
law interpolation between the sub-mm and X-ray wave bands [p7|], ignoring the observed break at optical frequencies; this is also 
justified by the observation that in flares the important optical-to-X-ray spectrum of low luminosity blazars seems to be flatter 



than the typical (aox) ~ 2 seen in the quiescent state [ [40| , |50|j51| ]. In the following, we adopt a = | and introduce the Doppler 
factor and magnetic field in canonical units, Vi = V/IO, and Bi = B/10 G. For the break energy we use ^ 10-4 eV @], 
and e ^ 1— lOkeV from X-ray observations of flares in Mkn421 and MknSOl [40 53], leading to 6ph ~ 10'^. The luminosity 



at the spectral break is not taken from observed fluxes at eb, because the emission at low energy is likely to be superposed by 
the emission from other jet regions not associated with the flare; rather, we use the observed, isotropized X-ray luminosity at 
£x ~ 1 keV of the flare, Cx — [lO'*^ erg/s]£x,45, and determine Ch from scaling with the assumed power law photon spectrum, 
£b = 'Cx,45(eb/ex)^^'^; note that feph'Cb ~ 3£x- Opacity requirements suggest that for low luminosity TeV blazars the external 
radiation does not dominate over the synchrotron component in the comoving frame of the flow [^]; if this would be the case, 
the thermal-like properties of the disk radiation would yield a neutrino spectra substantially different from our results. 

For the characteristic Lorentz factor in the comoving frame we find 7b ^ [3x10^^]!?, which is of the order the maximum 
proton Lorentz factors observed in cosmic rays. We will show below that acceleration of protons in blazars cannot reach higher 
Lorentz factors, and confine ourselves to the case 7p < 7b. In order to ignore the upper limit in the photon spectrum, we have to 
require 7p > lO^-'^Vi. 

To estimate the relevant time scale for the transient emission, T, we use the X-ray variability, where we have to consider 
electron cooling times of order [30 s] B^^^^V^^ in the observer frame, which for standard parameters are much shorter than 
the observed rise or decay times of the flare. This suggests that one can use Eq. (|l3|), which explains the generally longer decay 
times by Lorentz boosting effects in transrelativistic internal shocks, rather than by slow cooling. We note that this is in contrast 
to the usual interpretation of T as the doubling time of the flare, which assumes much longer cooling times expected in the 
weak magnetic fields required by purely leptonic emission models to explain > %. The latter explanation may also apply to 
hadronic models, if second order Fermi acceleration plays a role, so that in principle T ~ — % can only be considered as an 
upper limit. The typical time scales of blazar flares are then T = [10"* s]74 with 0.1 < 74 < 10, corresponding to a comoving 
linear size of the emission region ~ [3x 10^^ cm] T^Vi. The fact that we probably have to deal with transrelativistic shocks 
also suggests 9p sa 1, which we assume in numerical estimates; in general, however, we keep 9p as a free parameter. 

In a free relativistic jet we have iad ~ Tci, i.e., the emission region expands with the velocity of light, but /3cx < 1 is 
also possible if the jets are confined. Heuristically, we can express the value of P^x by the opening angle of the jet, (3cx ~ 
min(l, SjctFjet), and assume 8jct ^ (0viow)bi ^ 0-1' which means that Fjot ^ 10 corresponds to f3cx ^ 1- Hence, existing 
observations cannot decide whether magnetic confinement applies or not, but /3cx can be expected to be close to 1. On the 



other hand standard jet theory assumes that B (x R [54|, i.e., a — 1, but a — 2 may apply if reconnection isotropizes the 



magnetic field as assumed in Gamma-Ray Bursts [p5[]. Since (3cx and a always appear as a product in the equation, a/3cx ~ 1 
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FIG. 1. Dominant cooling processes and neutrino spectral shapes for AGN jets — Left: Parameter space, with the "star-point", denoting 
equal cooling time scales at the maximum energy, indicated. The shaded regions correspond to the dominant cooling process at the maximum 
proton energy: (A) Larmor limit or adiabatic cooling; (B) photohadronic cooling; (C) and (D) synchrotron cooling, where (D) marks the region 
where photohadronic cooling dominates for a part of the energy spectrum. Also shown are the positions of three observed AGN flares: (1) 
Mkn 421, April 26, 1995 (T4 = 10, £x,45 = 0.5) [|^; (2) Mkn 421, May 7, 1996 (T4 = 0.1, £x,45 = 0.9) |||]; (3) Mkn 501, April 16, 1997 
(74 = 3, £x,45 ~ 2.0) ||4c|]. Central positions assume ub = itph, black triangles correspond to ub ~ lOOtiph, diagonal errors indicate the 
range of possible Doppler factors (see text). Numbers in diamonds assiciate data points to the corresponding delimiting lines of muon cooling 
(black) and pion cooling (white); secondary particle cooling is relevant in the parameter space region above these lines. Right: Schematic 
representation of the shapes of neutrino spectra (time integrated power per logarithmic interval of energy). In Lv(Ev) vs. \uEv, corresponding 



to regions (A)-(D). Break energies due to changes of the dominant proton cooling process are indicated (cf. Sect. [I A), possible additional 
breaks due to secondary particle cooling are omitted for simplicity (cf. Sect. Ill C 2 1. 



is a reasonable assumption, and will be used in the following. The energy limit set by adiabatic cooling then corresponds to 
Pl ~ ^' so '^hat our canonical assumption ~ 1 is equivalent to the assumption that particles are accelerated up to their 
Larmor limit. Using these standard parameters, we find for the star-point of the parameter space 
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Fig.[i| shows the different regions of dominant cooling at the maximum energies, and their associated spectral shapes, where 
we have scaled Vi and Bi relative to the star-point values. Also shown are the positions of three observed AGN flares in 
the parameter space, for which we assumed the magnetic field to be in equipartition (a) with the radiation density, or (b) with 
the energy density of protons, Up = lOOiig, corresponding to standard assumptions in hadronic blazar models, and a range of 
possible Doppler factors 0.3 ^ Vi ^ 3. Comparison with the regions of dominance of muon and pion cooling corresponding 
to these flares shows that muon particle cooling plays a role for short flares in scenario (b), if I?i ^ 1; pion cooling is mostly 
unimportant for usual hadronic AGN models. In most cases, the neutrino energy is limited by Larmor radius constraints of the 
accelerated protons, consistent with earlier assumptions [|1 ll|16|]. 



3. Blazar neutrino maximum energies and fluxes 



Keeping the magnetic field, the Doppler factor and the proton-to-e lectro n energy ratio in blazars as free parameters, rather than 
adopting common assumptions, we can apply our discussion in Sect. Ill D| to obtain general upper limits on neutrino energies and 
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fluxes from this class of objects. We see from Fig. |l] that, for parameters typically observed in blazar flares, that Bsyn(fi) (T^) < 
^syn(p) (2^). but -Bsyn(7r) (2?) ^ ^syn(p) (2?) fpr all Dopplcr factors in the discussed range. Therefore, £^ < £^,p {T>) applies to 
muon neutrinos from pion decay, and £1, < £^_^{T>) to muon and electron neutrinos from muon decay. 

To find out the relevant range for the Doppler factors, we start with energetical considerations. The usual limit applied to the 
power of AGN jets is the Eddington luminosity of the putative Black Hole in the AGN, ijot ^ [10''^^ crg/s] A/bh,9, where A/bh,9 
is the mass of the Black Hole in units of lO^M©. Since we consider a beamed emitter, we have to set Q, = . Inserting in 
Eq. (^6|), assuming an energy dissipation efficiency of ^ 0.2 in the jet as expected for transrelativistic internal shocks, and 
equipartition of proton and magnetic field energy density, = 1, we obtain ^_B-y,max ^ [4x 10"^] A/bh,9'Cx 452^1, thus for the 
limiting Doppler factor allowing fj^^T = ty^-,^, 

Pi,, « 2 X M%T-^l^d-^l^ . (40) 

For the canonical range of assumed AGN Black Hole masses, 0.1 < Mbh.9 ^ 10, we can therefore assume that > T)* , for 
which case we find 

£vAT^i) ~ 4xl0i*^eV X D^^'^rll'^d''^^'^ < IxlO^^eV x M^^^^Q-^ (41a) 
LAT^i) ~ 0.3 X Dl^'Ti'/'el^' < 0.3 X Mli'^T.^'/'e'J' , (41b) 

where we assumed bp — 20. Adopting these parameters, neutrinos from muon decay are hmited by muon synchrotron cooling 
to an energy ^ [SxlO^^eV] M^^^^gXl^^e;^^'^ . On the other hand, we can also find the conditions under which neutrinos 
from muon decay, in particular electron neutrinos, can reach their highest energies and fluxes. We determine syn(n) ~ 
3 /:i(45r//^6'y^ leading to 

Pi,^ « 5 X Mli^gT^-^^^'e;'/'' . (42) 

Again we can confine the discussion to the case 2?^ > T)*, which gives 

~ IxlO^^eV X Dj/V//^6lp^/^ < 1x10^^ eV x Ml(^^g0-^ (43a) 

4^(2^1)- 2x10-^ X Dl^^T^/^el^^ < 5x10-2 X A/'/^g0p/^ (43b) 

The result that the upper energy limits for neutrinos from pion and muon decay are equal in their respective optimization is a 
consequence of the result that th e maxi mum energy of pion neutrinos for Doppler factors V > 2?, remains unchanged at the 
value Oy^TT, as discussed in Sect. QID 2. For the assumed range of AGN Black Hole masses and 6p > 1, we therefore obtain a 
strict upper limit of about 3 x 10^^ eV for neutrinos from AGN flares, which is independent of the flare time scale and any model 
assumptions. The inverse linear dependence on 9p, however, shows that this energy limit would strongly increase if we assume 
acceleration on time scales much shorter than the particles Larmor motion. 

If the TeV emission from blazars, which has generally a luminosity comparable to the X-ray emission, ought to be explained 
by hadronic emission, the corresponding neutrino luminosities would have to be of the same order, or ^ ^. Eqs. ( p3| ) show 
that this is incompatible with the conditions neutrinos from muon decay (electron neutrinos) need to reach their theoretical 
energy limit of ^10^^ eV. For muon neutrinos from pion decay, however, this seems to be possible for very large baryonic and 
magnetic energy densities, ^B-y ^ 10"^; this scenario would expect a difference between the electron and muon neutrino cutoff 
energies of more than one order of magnitude. It should be noted, however, that proton blazars can produce Cy{£u) ^ £x also 
for relatively moderate values of £,bi, if I? ~ P*; then, the hadronic radiative efficiency (neutrinos + gamma-rays) increases to 
<;50%, so that a comparable emission in neutrinos, high energy and low energy photons can be achieved fo r ~ £. b^ ^ 1- 
This more realistic scenario leads to maximum neutrino energies much below the upper limits stated in Eqs. { Wz) and (4-3a), so 
that flares from AGN jets would not be expected to emit considerable neutrino fluxes above a few times 10^^ eV. 

To get an estimate on event rates in current or planned VHE/UHE neutrino observatories, we consider the example of the May 

7, 1997 flare of Mkn 501, which lasted 3x 10* s: The total isotropized energy emitted in optical to X-ray photons of this flare is 

£ox ~ 2x10^° erg; taking the luminosity distance of wl60Mpc, and assuming Cy ^ Cox as suggested by TeV observations, 

this corresponds to a total energy in neutrinos of £u ^ lO^'^eV at earth of about 6x 10^ ergkm"^, which would produce 

'--^Sx 10^^ neutrino induced showers per km'^ air volume. The biggest air fluorescence detectors currently planned would cover 

about 3xl0^km'^ air, and the ground array of the Pierre Auger Observatory would correspond to about SxlO'^km'^, which 

2/3 

would clearly be not sufficient to detect the neutrino emission of a single AGN flare. Following an £u ' power spectrum, the 
total energy in neutrinos at £y ^ 10^^ eV would be about 6x10'^ ergkm^^, which would cause ^3x10-'' events in a 1 km"^ 
underwater/ice Cherenkov detector. Also here, even the biggest neutrino telescope currently considered would not be able to 
"see" single AGN flares. The best we can expect is therefore to collect diffuse fluxes corresponding to many AGN flares and 
determine the average properties of their neutrino spectra. Details of the time integrated emission spectra of AGN correlated 
transients, however, would still be important to determine reliable estimates for such diffuse fluxes. 



17 



B. Gamma-Ray Bursts 



Gamma-Ray Bursts (GRBs) are thought to be produced in highly relativistic outflows originating from a compact, explosive 
event over time scales of less than a second up to several minutes 1 57 ,58]]. Recent observations of ORB afterglows indicate 
that they are located at cosmological distances, which requires a characteristic luminosity of about ^lO'^^ ergs~^inder the 
assumption of isotropic emission. Possible scenarios which could release that much energy over the time scales observed are, 
e.g., the coalescence of a neutron star binary, or the collapse of a supermassive star The radiation observed from GRBs is 
expected to be due mainly to synchrotron or inverse Compton radiation from relativistic electrons accelerated at shock fronts 
occurring near the interface of the expanding relativistic shell (external shocks), or at shocks forming within the unsteady outflow 
itself (internal shocks) [|60[|. The same mechanisms would also accelerate protons, which could reach energies of the order of 
the highest energy cosmic rays, ~10^° eV [ |6l|j6^ ]. Due to the interaction with the dense photon field in the burst, these protons 
can produce efficiently VHE, and maybe also UHE neutrinos Jl4|]. Neutrinos of lower energy may also be produced by pp 
interactions between cold protons in the colliding ejecta [[l3|]. Obviously, both scenarios fall into the class of transient emission, 
and we can apply our results from Sect. |l| to examine in more detail the spectrum and maximum energy of the neutrinos 
from p7 interactions in Gamma-Ray Bursts. We note, however, that our discussion assumes that the physical parameters in the 
transient remain approximately constant over the emission time scale; it does therefore not apply to GRB afterglows, in which 
the parameters change drastically over very long time scales. 



1. Cosmological fireball models and internal shocks 

Many GRBs show intrinsic variability on time scales T ^ 1 ms— 1 s, while the total burst durations are typically T^rtB ^ 
0.1— 100 s [^. This implies that their energy is released in a volume of the typical dimension of compact or stellar objects, 
Rq ^ 10^— 10^° cm. For a total energy of >10^^erg, this leads to a local photon density of >10^^ergcm~^ with photon 
energies ^rrieC^. These "fireballs" (and even much weaker ones as well) would be optically thick to 77 pair creation, and 
for small baryonic loads. A/bar ^ 10~^Mq, the expansion leads to a conversion of almost all the radiation energy into bulk 
kinetic energy of motion, accelerating to a limiting Lorentz factor of F ~ ?7 = CcRB/A/baiC^ ^ 1, before photons can escape 
[|63|]. Hence, a dissipation mechanism reconverting the bulk kinetic motion into radiation is required after the flow becomes 
optically thin: this can be achieved by electron acceleration at shock waves occurring when the ejecta run into external matter 
[|63|]. Moreover, internal shocks can form in the ejected wind if the outflow is non-steady, i.e., 77 varies significantly over time 
scales ^Tgrb, which can lead to faster shells catching up with slower shells [04], similarly to what was discussed for AGN 



in Sect. Ill A . Such internal shocks can dissipate the kinetic energy with an efficiency comparable to external shocks. While 
external shocks are expected to produce a relatively featureless outburst over time scales comparable to the total burst duration, 
T^RB, internal shocks could be associated with the rapid variability within the burst on time scales T ^ Tgrb- In both external 
and internal shocks a substantial fraction of the gamma radiation is produced by synchrotron cooling of the shock-accelerated 



electrons. At the dissipation radius where internal shocks occur [241, 



rd ^ cTt]^ , (44) 

i.e., where the radiated gamma-rays are produced, the comoving magnetic field energy density can be parametrized through 
Ub — S,B-yUph, SO that 

B ^ lO^^Gx (ei3^6ph£5i)'/'VN"' . (45) 

where Tq — T/1 s and £51 — £/10^^ ergs^^ are the normalized GRB variability time scale and isotropic luminosity at the 
break energy in the observer frame, respectively, and feph is the bolometric correction factor correlating this specific luminosity 
to the total gamma-ray luminosity of the burst. The value of 6ph depends essentially on the high energy cutoff in the photon 
emission, which cannot be inferred from current data; for internal shocks one expects e < 100 McV in the comoving frame, 
corresponding to a canonical value feph ~ 10. In the internal shock model, constraints on the bulk Lorentz factor can be inferred 
from the requirement that the dissipation radius is larger than the radius of the photosphere, Tph ~ [10^^ cm] £5177"^, below 

which the wind is optically thick, and the radius of the external termination shock, r^sh ~ [10^^ cm] £g(^7]j^^^gr/~^/'^, which 
requires the bulk Lorentz factor of the ejecta, F ^ 77, to be in the range 30 < F < 10^, leading to magnetic fields in the range 
10^ G Jc; B Jc; 1 G, assuming equipartition between magnetic field and photons, standard values for Ch and foph, and typical 
time scales of 0.1 s for short-term variability, and 30 s for the total duration of a featureless burst. The latter values, F ^ 10'^ 
and i? ~ 1 G, would imply that the internal shocks occur on similar time scales and physical conditions as the external shock, 
making both scenarios virtually identical with respect to the discussion of transients. Typical parameters for assumed internal 
shocks with T <C T^rb are F ~ 300 and B > lO'^ G, which are also required by models predicting the acceleration of UHECR 
protons in this scenario [pTlJ 
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We still need to relate the parameters R, and V, defined in Sect. [II B for a causally connected emission region moving at 
some angle 9vicw to the line of sight, to the parameters describing the expanding flow in a GRB. In the comoving frame of the 
wind (at an arbitrary point), the apparent thickness of the wind zone extending to a radius r is r/T. Similarly, the transversal 
extent of causally connected regions in the comoving frame of the flow is r/T, because regions farther apart move away from 
each other with velocities larger than c {r/T may thus be interpreted as the "Hubble radius" of the expanding emission region, cf. 
[]63|]). If the energy dissipation takes place at a radius rd where the bulk Lorentz factor is saturated, F ^ 77, the comoving "linear 
size" of the emission region can thus be written as i? ~ rd/77. On the other hand, the isotropic luminosity of the burst in the 
observer frame is related to the photon energy density in the comoving frame of the wind by £ 47rr^cr^7iph, which is identical 
to Eq. ( p^ in conjunction with the linear boosting formula, C = LV^, if we use the parameters R — r^/T, = 1 and V — T. 
For T ^ T] ^ const, all parameters are therefore correctly derived if we treat the GRB as emission from a spherical region of 
radius R ^ r^/rj boosted with a Doppler factor T) ^ rj. Moreover, if the gamma-ray emission originates from synchrotron 
radiation of electrons whose cooling time scale is much shorter than the crossing time [p4|], we can also apply R « cTV, or 
Pt ~ 1; since we can apply our discussion to both internal spikes in strongly variable bursts and to featureless bursts in total, we 
allow the value of the normalized time scale in the broad range 10^'' < Tg < 100. According to our definition of the radius R as 
the "Hubble radius" of the expanding emission region, we obviously also have /3ox = 1, and the effective space fraction available 
for the gyration of protons is essentially constrained by adiabatic losses, ~ (7ra6'F)~^, where both a = 1 and a = 2 have 
to be considered possible values, depending on whether the magnetic field is largely transversal or isotropized by reconnection 



2. GRB neutrino spectrum and maximum energy 

The observed electromagnetic spectrum of a GRB can be approximately described as a broken power law, with a break energy 
£b,GRB ~ 300 keV in the observer frame. The photon number spectrum is then given by A^(e) « for e > eb,GRB (thus 
a = 2), and for e < eb,GRB 

it is N{e) 05 e"^/^, with 

£b,GRB — EbjGRB^ ^ ■ This yields a proton break Lorentz factor in the range 
^ 7b 5x 10^, corresponding to the range of possible Doppler factors given above. When we consider neutrinos of energy 
£1, ^ 100 TeV, we require proton Lorentz factors of -fp > 3x10^/1). For V ~ 100, we therefore have 7^ > 7b; for simplicity, 
we restrict our considerations to the case 7^ ^ 7b, noting that this might be only marginally correct for the lower energy bound 
of the VHE neutrino regime. In a e^^^^ low energy spectrum, we simply have A^ph/^ph.b = 4 — 3(7b/7p)^/'^ for 7p > 7b; 
for simplicity, we use iVph/A^ph.b ~ 3 for all 7^. Introducing normalized quantities also for the Doppler factor, V = 100 1?2, 
the magnetic field, B = [lO"^ G] B3, and the bolometric correction factor, &ph — IO61, the coordinates of the star-point of the 
parameter space, where all proton cooling processes have equal time scales at the maximum proton energy, are then 

0.9 X Cli%-'^'a-'^^ (46a) 
4.5 X £^^^%'^^^dl^^a'^/^ (46b) 
2x10-2 X Cli'T-'^'el/'a'H^' . (46c) 

In the following we express the magnetic field by its equipartition parameter. Figure || shows the GRB parameter space and 
the separate regions of dominance of the various cooling processes limiting the neutrino energy, and the corresponding neutrino 
spectral shapes for a set of possible parameters, including both millisecond flares in internal shocks and featureless GRBs; unlike 
in the AGN case, we do not use parameters scaled relative to the star-point, and put more emphasis on secondary particle cooling, 
which we consider separately for pions and muons. We see that in all cases, the energy of neutrinos from both pion and muon 
decay is limited by secondary particle cooling, where synchrotron cooling of pions and muons plays the most important role. 
Only in a limited part of the parameter space does an additional break due to adiabatic cooling of muons appear in the spectrum, 
while adiabatic cooling is generally unimportant for pions. At the star-point, the maximum neutrino energy from Gamma-Ray 
Bursts is therefore considerably below £*, so that we omit its value here to avoid confusion. 

Because of the dominance of secondary particle cooling, Eqs. ( pT| ) must be used to determine the neutrino energy limit at a 
given Doppler factor, £^{1)), which is reached along a line = S/g^ syn(*)(^/^*)^ ™ parameter space, with 

eB,,sy„w = 2.3x10-4 y. d,{'e,T-''\-^'%^^ (47a) 
^s^.sy„(p) - 1.7x10-5 X . (47b) 

Recent observations of GRB afterglows at large redshifts require a total isotropic energy emitted in photons of 6ph'Cb'?GRB ^ 
2953 ej-g Since this value comes very close to the gravitational energy released in the collapse of a compact object (e.g., 
a neutron star), it is most likely that the energy of a GRB is not emitted isotropically. Moreover, efficient hadronic emission 
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FIG. 2. Dominant cooling processes and spectral shapes for neutrino production in Gamma-Ray Bursts — Upper left and lower panel: 
Parameter space for three different parameter sets, all assuming Op — bi — 1, and a = 2: (A) canonical case, £51 — To = 1; (B) short intrinsic 
flares, T = 10 ms, and £51 = 1; (C) extreme case for bright afterglow burst, £ — 3x 10^^ ergs~^ and T = 30 s (which implies an isotropic 
bolometric photon energy bphChT ~ 10^^ erg, requiring <^ 1). The central line divides regions where photohadronic cooling dominates 
adiabatic cooling (left from the line), and vice versa (right from the line), the star-point of equal proton cooling times at the maximum energy is 
indicated. Shaded regions correspond to spectral shapes produced by the subsequent change of dominant proton and secondary particle cooling 
processes: (1) adiabatic/photohadronic coohng dominant up to 8^; (2) £'y,ad(*) < £v', (3) f^,ad{*) < £v.Byn{*) < £-v', (4) f,/,syn(*) < £-v\ (5) 
£}^,syu{*) < £syn{p) < fiy. In cach figurc, the upper part corresponds to neutrinos from pion decay and the lower part for neutrinos from muon 
decay. Also indicated are the regions allowing UHE cosmic ray production up to f p,20 (black chain lines), and the region corresponding to the 
neutrino fluxes predicted in Ref. |[l4[], for = 1 (white hatched line and arrows). Upper right: Schematic representation of corresponding 
spectral shapes, log £^ (£v) vs. log corresponding to regions (l)-(5). The lower break indicates the spectral change at about 100 TeV due 
to the change of the photon target spectrum at 7p ~ 7b (see Ref. [O]), which is not discussed in this paper. 
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of GRB has to assume that comparable amounts of energy are present also in other channels, like magnetic fields or protons. 
While collimation into jets is a distinct possibility, the evidence for it is not as obvious as in AGN; in the following we assume 
that the GRB energy is emitted into a "firecone" of solid angle Att ^ fl ^ "D^^, and we introduce the normalized solid angle 
ri47r = il/4TT ^ 1. For GRB with a strongly variable light curve, the transients considered here are single, isolated radiation 
spikes rather than the whole burst; if we assume, however, that the total energy is equally distributed over the individual flares 
and that there are no extended gaps in the light curve, we have £t /bphCh"^ ~ ^grb / bphChTc^B, which again allows a common 
description of both subflares within GRB and featureless bursts. For simplicity, we assume in the following parameters which 
allow relativistic protons and magnetic fields to dominate the energy, £^pB ~ 1 and ^57, max ~ ^4^^5i ^ neglect the 

dependency on the weakly varying factors a and 61, leading to 

1)2,. « 5 X n^y'T'^/'e;^^' (48a) 
I>2,M « 8 X n~y%~'/'9;^^' , (48b) 
which means that we only need to discuss the case T>^ > V* . We then obtain for the neutrino energy limit 

L,^{Vi) w IxlO^^cV X dI^%^/^0-^^^ < IxlO^^eV x n^^^^To^^e-^^^ (49a) 
« 4xl0i^eV X dI^\'/^9-^^^ < IxlQi^eV x n~y^"Ty''9-^^' . (49b) 

It therefore seems that photohadronic neutrinos from GRBs can reach energies up to 10^^ eV and above. However, the corre- 
sponding normalized neutrino luminosities at the neutrino energy limit are 

L^CD2) ~ 2x10-6 X Dae/p < 1x10-5 X 17^/^ V^/^6lp/^ (50a) 
~ 2x10-7 xI?20F < 1x10-6 X rj^/^V^/^^y^ , (50b) 

which shows that even in the limiting case T) = T) only very small fractions of the photon energy can be emitted in neutrinos 
at the limiting energy. Since we assumed ^ 1 and ^s^^max Jc^ 1 in our calculations, the low values of are clearly due 
to a very low efficiency along the line = 6_B7,syn(*)(2')- Assuming a neutrino conversion efficiency of C,^ ^ 0.2, S^B-y ^ 1 
and Up ^ 5e^((iA'ph/(ie)b, which is consistent with our equipartition assumption, ^ S^bi ~ 1, for feph ~ 6p « 5, viz., 

^ 0.04, Waxman and Bahcall [|l^ derived a GRB-related neutrino event rate at > 10^'* eV of about 10—100 per year in a 
km^ detector This event rate would be below the expected background, but still statistically significant because of the time and 
directional correlation to the bursts. It is clear from the above that if we try to maximize the neutrino energy in the UHE range, 
following Eq. (p9|), this estimate would have to be reduced by more than 4 orders of magnitude, leading to insignificant event 
rates at these highest energies even if the usually larger detector volume of UHE experiments is taken into account. 

We can also turn the question around and ask which is the region of parameter space where neutrino fluxes of the magnitude 
predicted by Waxman and Bahcall [14| are expected, i.e., > 0.04, and then derive the maximum neutrino energy for these 



parameters. The assumptions of Waxman and Bahcall for a single burst (or sub-flare) can be essentially put in the form Ci, ^ 
0.2ChT, which can be rewritten in the form of Eq. ( ^ as bphCi^^B-y^pB ^ 0.2bp. To optimize the efficiency, we have to evaluate 
Ci/ at an energy where it is not yet diminished by synchrotron cooling of either protons or secondary particles; since T) > V* 



this means that we have to use Eq. (32b) for a = 1. For hp ^ 6ph and ^pB ^ 1 this leads to the condition S^b-i ^ {'D/'D*Y for 
T) > V*; for V < V*, the efficiency is constant since photohadronic cooling dominates. Figure^ shows that this region of the 
parameter space is entirely enclosed by the regions of pion, muon and proton synchrotron cooling dominance, which means that 
Su,syn{*) has to be used as the maximum neutrino energy, according to our assumption of adiabatic coohng dominance. Inserting 
into Eq. (|o|) we obtain 

< IxlO^^eY xV^Ta^^ < 6x10^^ eY x Q^^^ (51a) 

< 5xl0i5eV X P^T-^^/^ < 4x 10^^ eV x f^-J^^ , (51b) 

where the second limit was derived using S^B-y < ^41^^51 ^"d il47r <C 1. This means that Gamma-Ray Bursts can still produce 
UHE neutrinos with considerable fluxes, but probably not above 10^^ eV. 

While this paper was in preparation, it has been proposed by Vietri [^sjj that neutrinos of more than 10^^ eV can be produced 
in external shocks of GRB and in GRB afterglows, with fluxes observable in very large scale air shower experiments. Concerning 
the external shock of the main burst, this prediction is clearly in conflict with our upper limit stated in Eqs. ( ^0| ) and (|5l|). In the 
afterglow, the larger time scales and lower photon energies make the situation rather comparable to blazar jets, and Eqs. ( |5l] ) do 
not apply; on the basis of the considerations presented here, we cannot rule out the possibility of producing neutrinos of such 
extreme energies in afterglows. It has been recently shown that such large UHE neutrino fluxes would also be in conflict with the 
assumption that (a) the cosmic ray production spectrum is cx 7-^, as assumed also here, and (b) that the locally observed cosmic 
ray energy density above 10^^ eV is homogeneous throughout the universe and does not evolve with cosmological redshift [tea]. 
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3. The relation between cosmic ray and neutrino maximum energies 



GRBs have been proposed as possible sources for the highest energy cosmic rays [ piy62| ], which are observed up to £p ^ 
3 X 10^" eV [|9|] . The maximum proton energy in the observer frame is £p = nipC^jpI); at the star-point of the parameter space, 

we find £* ^ [Ix lO'^^ eV] C\{'^9p '^^'^ a~'^l'^ . The highest proton energy for a given Doppler factor is obviously achieved at the 
line of equal time scales for proton synchrotron colling and adiabatic losses (the same condition which determines the maximum 
neutrino energy if secondary particle cooling plays no role). This is the border line of region (5) in the parameter space shown 
in Fig. H so we find Ep = £*{V/V*Y/^ for V>V\ leading to 

£p{V2) ~ l-x.lQ^'^cW xVl'^T^'^e^'^'^a-^/'^ . (52) 

Since ^b-^ oc P^''/'^ along this line, the required magnetic field equipartition parameter rises fast; as a function of the proton 
energy in the observers frame, 4,20^ in units of 10^° eV, we can formulate the minimum requirements in V and as 

:Puhe(4,2o) « 90 X £l%%-^'^el'^a^'^ (53a) 
eB7,uHE(4,2o) « 0.01 X £l%T-^^^ela^'^C^^b^^ . (53b) 

In correspondence with the result of Waxman [61 1, we find that bulk Lorentz factors T) > 300 and magnetic fields close to 



equipartition with the radiation (^s-y ^ 1) are sufficient to produce the highest energy cosmic rays, 4,20 ^ 3, provided that 
protons are accelerated on their Larmor time scale. Assuming minimal conditions for the production of UHE cosmic rays, i.e., 
T) — 2?uHE(fp,2o) and ^b-^ = £,B-y,v-a-E{£p,2o), we can derive the corresponding maximum neutrino energy as a function of £paa 
as 



f 

f 



£p=£p(r>) 

£p=ip{V) 



SxlOi^eV X 4/>y Vo^/^a-i/^ (54a) 
3xl0i6cV X 4/>y (54b) 



where we again assume £^ i, = £^^syn{*)- We see that, although secondary particle cooling limits the neutrino energies to values 
much below the canonical 5^ = 0.05£p assumption, the maximum neutrino energy still has a tendency to increase with £p. 



However, Eqs. (54) assume just minimal conditions for the production of UHE cosmic rays, while the region of the parameter 
space corresponding to £p > [lO^*' eV] £p,2o, as shown in Fig.^, allows neutrino break energies principally somewhat above or 
below the value stated in Eqs. ( |54| ) for given £p,2o- 

We note, however, that according to the standard assumptions of shock acceleration, neutrino production and cosmic ray 
ejection from magnetically confined sources are physically connected processes. Although it is a distinct possibility that cosmic 
ray acceleration to £p > lO^*' cV on the one hand, and efficient VHE neutrino production requiring only cosmic rays of lower 
energy on the other, may happen at different radii in the expanding shell [^], the ejection of UHE cosmic rays is non-trivial in this 
scenario, since the particles are advected downstream, and thus accumulate inside the expanding shell and remain magnetically 
confined. As the shell continues to expand and the comoving magnetic field decreases, the cosmic rays would lose most of 
their energy by adiabatic cooling, before they eventually escape when the shock slows down. The easiest way to circumvent 
this problem, and to eject cosmic rays with the high energies they receive at the shock, is to convert them into neutrons in, e.g., 
P7 — + n7r+ reactions, which decouples them instantaneously from the shell as long as their reconversion probability is low 



|]64|], which is the case for T> > T>* (cf. Sect. [I A). This means that a low neutrino production efficiency corresponds to a low 
cosmic ray ejection efficiency for the standard scenario of shock acceleration, assuming magnetic confinement of all accelerated 
charged particles. This physical connection of both processes makes the hypothesis that Gamma-Ray Bursts are the sources of 
the highest energy cosmic rays testable by neutrino VHE and UHE observations, both with respect to the flux and the maximum 
energy of the putative ORB correlated neutrino spectrum. It is also obvious from Fig.^ that this joint scenario, i.e., 4. 20 > 3 
and £^ > 0.04, also requires extreme values of the GRB parameters V and £,bj- 



V. CONCLUSIONS 

We have presented a detailed investigation of the production processes of very energetic {>1Q^^ eV) photohadronic neutrinos 
in relativistically boosted astrophysical sources. Using the constraints set by the source variability, and assuming that the 
acceleration process for protons is of the Fermi type, we derived limits on the maximum energy and the position of possible 
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breaks in the neutrino spectrum. Comparing the effects of various cooling processes for both protons and secondary particles 
in the hadronic cascade leading to neutrino production (i.e., pions and muons), we find a general upper limit on the neutrino 
maximum energy, which does only depend on the Doppler factor of the emission region relative to the observer. Energetic 
constraints allow one to turn this into a general upper limit, which is only dependent on observational parameters of the transient, 
but not on any model dependent parameters. This is the major result of this paper. In some cases, and assuming that the energy in 
protons, magnetic field and photons are near equipartition, stricter limits can be imposed when considering both neutrino energy 
and expected flux. We apply this general result to two classes of proposed cosmic neutrino sources: hadronic models of Doppler 
beamed jets from Active Galactic Nuclei (AGN), also called blazars, which are known to emit most of their energetic radiation 
in short, distinct flares, and Gamma-Ray Bursts (GRBs). 

For blazar flares, we confirm that under the most common assumptions the neutrino energy is limited to ^10^* eV by Larmor 
radius (or adiabatic) constraints of the accelerated protons. For short (<1 hrs) flares, however, the maximum energy of neutrinos 
from muon decay may be additionally suppressed by muon synchrotron cooling. Assuming that the AGN is fueled by Eddington 
limited accretion on a supermassive Black Hole with A/bh < we show that neutrinos from AGN flares cannot exceed 

energies of ~10^^ eV, independently of the time scale of the flare. For electron neutrinos, which result only from muon decay, 
we can show that fluxes of the magnitude usually assumed in the literature (i.e., similar to the X-ray luminosity of the source) 
can only be attained if the energy extends only up to about 10^^ eV. Unless vacuum neutrino oscillations occur in nature, this 
has important implications for the neutrino event rates expected in the Pierre Auger Observatory, where electron neutrinos are 
proposed to be most easily detected because of the distinct properties of their induced air showers. 

For GRBs we find that the synchrotron cooling of pions and muons limits the maximum neutrino energy over most of the 
allowed region of parameter space. We show that, although neutrinos from GRBs are in principle able to exceed 10^^ eV, in 
particular for acceleration over long time scales in or near external shocks, this possibility would imply extremely low efficiencies 
and thus insignificant neutrino fluxes. If we require that neutrinos are also produced with fluxes similar to the gamma-ray flux, 
and applying the usual energetic constraints for near-isotropic GRB sources, we find an upper limit on the neutrino energy of 
<10^^eV for muon neutrinos (from pion decay), and <10^^eV for electron neutrinos. This limit can only be increased to 
UHE (>10^^eV) neutrino energies if strongly coUimated outflows are assumed. We also show that the conditions for GRBs 
to accelerate protons up to the highest energies observed in cosmic ray air shower experiments (^3xlO^°eV) coincide with 
the conditions for efficient neutrino production, and expect a flat neutrino power spectrum extending up to a break energy in 
the range of 10^^—10^^ eV. If neutrino production and cosmic ray ejection from GRB are connected processes, as implied by 
the standard assumption of magnetic confinement of shock accelerated particles in an expanding shell, this would make the 
hypothesis that Gamma-Ray Bursts are the sources of the observed ultra high energy cosmic ray spectrum testable with neutrino 
observatories. 

As a corollary of our investigation of the relevance of secondary particle cooling to hadronic cascades in common models 
of astrophysical transients, we have also checked the effect of this process on other proposed, non-transient sources for cosmic 
neutrinos. In particular the predicted diffuse neutrino background from AGN cores, frequently used for event estimates in high 
energy neutrino detectors, was previously derived disregarding secondary particle cooling. Here we obtain for energies above 
10^^ eV a strongly reduced contribution and a lower cutoff in the electron neutrino component, and a reduction of about a factor 
3 for the expected diffuse muon neutrino flux (see Appendix^). 

In cosmic sources where the neutrino energy is limited by secondary particle cooling, which is clearly predicted for Gamma- 
Ray Bursts and is possible in some blazar flares, the expected difference in the cutoff energy of electron and muon neutrinos 
could also serve as a laboratory to test the existence of neutrino vacuum oscillations in nature at very high energies — a detected 
change in the neutrino composition near the cutoff energy could rule out this possibility. Such a measurement would require a 
large detector sensitive in the range 10^^ — 10^^ eV, capable of detecting both electron and muon neutrinos and able to distinguish 
between flavors. 
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APPENDIX A: NOTATIONAL CONVENTIONS 

Below we explain some general conventions we use throughout the paper Table | lists some generally used symbols; it does 
not contain symbols which are used only in the section where they are defined. 
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TABLE I. Index of frequently used symbols. 



symbol 



meanmg 



definition/relations 



introduced in 



, r,, 7. 



m, 

A''ph, A*'ph,b, Eb, a 

RF 
^th 

7b, 7b 
tb,w. Ha 

B, ub,, 

tp 

fpi^ /syn, /ad, /ma 

7p, s 

r, /3r, Ovicw, 2? 

T, Lb, bph 

Lcr, ^nj, Trad 

Pt, a^t 

^acc , ^F 
a, /3ex, 

^;>,syn(p), ^u^py 

^;y,ad(*), ^;y,syn(*) 
EiRF 

it/,. (2?), -Bayn(.) 

Ct, C-B7,max 



particle mass, classical radius, and Lorentz factor [•=p,7r,n] 
unstable particle lifetime [•=ii,7r,n] 

Lorentz factor of secondary particles at production, decay'' 

total photon density, density and power law index above eb 

threshold photon energy in proton RF for 7r-production 

characteristic Lorentz factor for power law approximation 

pion production cooling time for 7^ = 7b, effective inelasticity 

weighted cross section for power law photon spectrum 

magnetic field, particle cyclotron frequency 

particle synchrotron [•=7r,n,p] and adiabatic cooling time^ 

total and specific photohadronic cooling times'* 

total proton cooling time scale 

rate of proton cooling relative to pion production 

efficiency for neutrino production, from specific decay'' 

total injected proton energy density, bolometric correction factor 

maximum proton Lorentz factor, power law index 

neutrino emission spectrum, cutoff energy, local spectral index 

bulk Lorentz factor, velocity in units of c, viewing angle and 

Doppler factor of the emission region in observers frame 

observed duration of transient, luminosity at e = eb 

transient crossing, proton injection and radiative time scale 

linear size of transient emitter, |i — in line of sight, ±= projected 

geometric correction factors 

proton Larmor radius and time 

acceleration time scale, normalized to Larmor time 

magnetic field decay parameter, expansion velocity of transient 

Larmor limit for neutrino energy, correction factor 

neutrino cutoff energy limited by proton cooling'' 

characteristic frequency for secondary particle decay'' 

critical neutrino energies for secondary particle cooling"'' 

neutrino energy in decay frame'' 

dimensionless characteristic parameters of transient 

magnetic, electron and photon energy density in emission region 

energy equipartition parameters 

"star-point" parameters'^, all specific proton cooling times —t^cc 
neutrino energy limit and corresponding B for given V [•=p,*]'' 
relative neutrino luminosity for £1, — £u.» 
total energy budget of transient, corresponding maximum 



r, = e 
r, = 7.r; 

7M=<7r«7:;^<7r~7p 

dN„h = e"" de, e > eb 
Eq.(g) 

7b = eSf/2eb,7b =7b/2? 

Eqs. (|), (|l|), (g) 
Wb,. = eB/m,c 
Eq.(|), tad = 2\B/B\ 

t"^ = t"' + t~' 

Eq.C) 
Eq.(|) 
Eq.( 



dN„ 



7p 



' djp, 7p ^ 7p 



Eqs.(§) 

/3r = Vi-r-^ 

D = [r(l-/3r COS Ovicw)] 
T = TV, Zlb = ibO^ Eq. (0) 
Eq.© 

= cTcr 

Eq. ^ 

rL = Ep/eB, ti^ = 2-Kri^/c 

ia.cc — Opti^ 

B oc R-", /3cx = R/c 
Eq.(ra,pL<min(|, ^g^l^J 
Eqs.^, ^ 

Eqs. (0) 

Eu,* = 7*jB^:* 
Eqs. (0), (0) 

Ub — B'^/StV, Ue ~ Uph 

Eq. Qi ^ps ^Up/uB 
Eqs.(|3|) 
Eqs.(Z7)-(0) 
Eqs. ( p), (H) 



general 
Sec 



Sec 
Sec 
Sec 
Sec 
Sec 

Sec 
Sec 
Sec 
Sec 
Sec 
Sec 
Sec 
Sec 
Sec 
Sec 

Sec 
Sec 
Sec 
Sec 
Sec 
Sec 
Sec 
Sec 
Sec 
Sec 
Sec 
Sec 
Sec 
Sec 
Sec 
Sec 
Sec 
Sec 
Sec 



OA 
[I A 



[I A 



[LA 
[LA 
[I A 



[I A 



[lA 
[lA 
[lA 
[LA 
[I A 
OB 
[IB 
[IB 

i 



[II A 1 

2l 



[II A 1 



[II B 
[II B 



[II B 


I 


[II B 


I 



[II B 



[lie 



[lie 



[II q 



[II Dl 



[II Dl 



[II D I 



[II D 2 



[II D 3 



[II D 3 



"Cooling processes are: adiabatic losses (ad), synchrotron radiation (syn), photohadronic interactions (p7), (charged) pion production (-n-.Tr*), 
Bethe-Heitler pair production (bh), neutron escape (csc,n). 

''The subscript "*" denotes secondary particles in hadronic cascade here and throughout the paper [*=7r±,M] 

"^The superscript "*" generally denotes quantities taken at the star-point of the parameter space, i.e., T) = D* , B — B*, and 7p — "/p. 



a. Units and normalized quantities: We use cgs units, except for particle energies which are given in standard muhiples 
of electronvolts (eV, MeV, GeV, TeV), and particle interaction cross sections measured in microbarn ( /ibarn). In numerical 
calculations we use quantities normalized to common powers of their standard unit (stu), X„ = X/10"stti (e.g., L51 = 

erg s^^). Dimensionless quantities may be normalized in common powers as well. This convention is used consequently 
in Sect. |^ which means that numerical subscripts always denote normalization powers. 

b. Reference frames: Three reference frames are used in the paper: the observers frame, the comoving frame of the emission 
region (relativistic flow), and the rest frame of an interacting massive particle. Quantities are by default given in the comoving 
frame; quantities in the observers frame are denoted by calligraphic letters (e.g., £, T, C, e) of the same kind as corresponding 
quantities in the comoving frame (E, T, L, e). Quantities in the particle rest frame are denoted with a superscript RF. 

c. Luminosity convention: Luminosities quoted in the paper always mean the isotropic radiated power at specific energy 
(frequency), e.g.,, L — e^{dN-p\^/ de) or L^, — El{dN^ / dE^). If we refer to bolometric luminosities, we do this by explicitly 
multiplying with a bolometric correction factor, e.g.Lboi = ^^ph^b (see Tab. |). The energy output of a transient over its time 
scale T is denoted the time integrated luminosity, L. 
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APPENDIX B: PHOTOHADRONIC INTERACTIONS 



The major photohadronic interaction channels of protons are single pion production with and without isospin flip, 737 ~+ nir* 
and p7 pir", followed by several two-pion production channels, and multi-pion production which dominates at very high 
interaction energies. Secondary neutrons can contribute negative pions from — > pn^ reactions, which are otherwise only 
produced in two-pion and multi-pion channels. The subsequent decay of the pions leads to neutrino production by 

7r± ^ f^^M'^f.) (Bla) 
e^v^Ve{v^Vf,) . (Bib) 

The charge of the initial pion is only relevant for the v : v ratio; it plays a role for the electron neutrino component at energies 
w 6x 10^^ eV, where the detectability of i/g is enhanced due to the W~ resonance in interaction with atmospheric electrons. 
Otherwise, neither underwater/ice, nor air shower experiments can distinguish between neutrinos and anti-neutrinos, so that 
we can disregard the sign of the pion charge. The average energy of the neutrinos is determined by decay kinematics: it can be 
written as (Ey) — E^^^j-^, where E^^^ is the energy of the neutrino in the rest frame of the decaying particle moving with Lorentz 
factor For pion decay, E^''^ = 30 MeV [^, while in the three particle decay of the muon {E^^^) « impC^ = 35 MeV; as 
an approximation, we may use {E^\) « jm.,^c^. 

At low interaction energies, e"'' ^ 340 MeV, the photohadronic cross section for both charged and neutral pion production 
is dominated by the Z\(1232) resonance, leading to single pion production with a 7r+ : tt" ratio of 1 : 2 following isospin 
symmetry, and £^7r ~ 0.2Ep from two-particle decay kinematics. These relations are often used as characteristic for pion 
production, and we call it the Z\-approximation [^|^] ; it is quite accurate for tt" production, in particular in steep photon spectra, 
and in thermal spectra with temperatures kT for proton Lorentz factors ■jp ^ [340 MeV] //cT, where higher resonances and 
other photohadronic interaction channels have little or no influence on the cross section. For charged pion production, however, 
other processes contribute significantly both above (e"'' > 400 MeV) and below (efjf < e"^"^ ^ 250 MeV) the dominance region 
of the Z\(1232) resonance, and thus enhance the contribution of charged pions relative to the Z\-approximation. These are in 
particular N* resonances at energies above the Zi(1232), but cross section data also require a contribution from non-resonant 
pion production, which give an almost constant background of about 100 /xbarn, extending from the immediate threshold (even 
before the Zi(1232) resonance becomes relevant) up to the highest energies [ p7|j6^ ]. In the neutrino sources considered in this 
paper, a proton spectrum extending up to a maximum energy Ep — rUpC^jp, interacts with an isotropic power law photon 
distribution, dA^ph c>c e"" de, a > 1, extending from a break energy Cb to some cutoff at e ^ eb, with a total number density 
^ph,b of photons with e > eb- Below the break energy, we assume that the photon number spectrum is flatter than e^^ 
everywhere. We can then define the inelasticity weighted effective cross section for pion production. 



Ha = 2 J dxx-" J2 ['^^ 



for a > 1 , (B2) 



where ef^ ~ 145 MeV is the threshold photon energy for pion production in the proton rest frame, ai is the cross section, and 
Ki — (AEp/ Ep)^ is the proton inelasticity of the reaction channel, evaluated at a photon energy e"-^ — xef^ in the proton RF 
(see for details). Introducing the characteristic Lorentz factor 7b = e{'j^/2eb of the protons, which is necessary to boost 
background photons at the break energy above the reaction threshold for pion production, we can then write the cooling time of 
protons with 7p = 7b as th,-K = [cAph.bHa]^^, where A'ph^b is the number density of photons with e > eb. As a function of jp 
the cooling time tp^^^ is then expressed by Eqs. (^, and we note that all the interaction physics, including the relative contributions 
of different resonances or other reaction channels, are absorbed in Hq and thus independent of the proton energy — this would not 
be the case in, e.g., thermal photon spectra, where our results are not applicable. As a numerical simplification, we also disregard 
the upper cutoff in the photon spectrum, which is justified if the spectrum for e < Cb is sufficiently steep — for spectra with 
a < 2, this approach is valid only in a limited range of Lorentz factors below 7b — and use Ha ~ H2/(a — 1) for 1.5 < a < 3, 
with H2 ~ 22 fihsLTii [|68|]. If jp 3> 7b, photohadronic interactions majorly happen at proton rest frame energies far above the 
threshold, where the cross section and efficiency is approximately constant; this justifies the approximation tp^j^ w [cA'phH2] 
used throughout the paper 

If no other cooling processes apply, the average number of pions produced per proton can be found as (N-^/Np) = lia/Ha ~ 



7, for 1.5 < a < 3, where Ha is defined as in Eq. ( |B2[ ) by replacing the proton inelasticity, Ki, by the pion multiplicity of the 
reaction channel. The average energy carried by each pion is then, independent of jp, given by {Et^) / Ep « i « mT^/nip, i.e., 
the Lorentz factor can be treated as conserved in the interaction, 7J' « 7^. Here, 7^' is the pion Lorentz factor at production, 
which has to be distinguished from 7^% the pion Lorentz factor at decay; for both pions and muons we consider the possibility 
that they lose energy during their lifetime, viz., 7^*^ < 7^'. Because of the small mass difference between pion and muon, we 
can also approximate 7^' w 7^^ Distinguishing between the charged pion and neutral pion multiplicity in the definition of Ha 
yields the charged pion fraction, tt* : tt" sa 2 : 1, almost independent of the power law index a [Bsl]. This result includes all 
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reaction channels, and is in contrast to the often used ratio 7r='= : tt" = 1 : 2, which is derived from the Z\-approximation. The 
discrepancy of a factor of 4 emphasizes the importance of charged pion production away from the Z\-resonance in power law 
target photon spectra, which is also relevant for the kinematics: it explains the difference of the usual assumption, {E-^) w \Ep, 
(which is strictly valid for tt° production, see above) and our result for charged pion production, {Et^) w \Ep. 

APPENDIX C: SENSITIVITY RANGE AND ENERGY DEPENDENCE OF NEUTRINO DETECTOR TECHNIQUES 

The low detection probability for neutrinos above the TeV range requires large detector volumes, which can be achieved, e.g., 
by the extension of classical water Cherenkov detectors to larger dimensions, such as the NESTOR or Lake Baikal (and the 
recently cancelled, pioneering DUMAND) experiments, and similar projects in the planning stage [^]. The same technique is 
also efficient in the deep antarctic ice, as shown impressively by the recent detection of the first neutrino events by the AMANDA 
experiment [^], and there is hope to extend this detector to an effective volume of 1 km"^ in the future [[7C| l. A cost-efficient 
way to further increase the detector volumes could be the detection of radio pulses [ [7l| ] or acoustic waves w% from neutrinos in 
water or ice, which however are limited to very high energies to obtain a reasonable signal-to-noise ratio [Ep!^ 

The neutrino event rate per logarithmic energy interval in a given detector, 0^,dct, can be written as the product of the 
neutrino number flux, Cu/£v, times the detection probability, Pu,dct- For deep underwater/ice Cherenkov detectors, which 
are most efficient to detect muons from ^ ^ conversions because of the large mean free path of muons in matter, the 
detection probability is essentially proportional to the ratio of the muon mean free path to the neutrino mean free path, yielding 
Pu,dct J5 ^v'^ for Ey ^ 10^^ eV. At energies above ~10^'' eV the effective solid angle covered by the experiment is reduced 
by earth shadowing effects [[73|], so that this technique becomes ineffective for ultra high energies; it is also less sensitive to 
electron neutrinos, because the short range of the electron in matter reduces the effective volume. Radio Cherenkov detectors 
are proposed to work best in deep ice. The detection probability is usually expressed as the effective volume of the detector, 
increasing as /V^dct for ^ 10^® eV, and is roughly constant above 10^^ cV [^T]]. The energy threshold for this 

technique is set by signal-to-noise constraints at ^5x 10^^ eV [|||, thus in the relevant sensitivity range the detection probability 
can be treated as approximately constant. No clear predictions exist about the efficiency of the acoustic method yet, which is 
probably most interesting at ultra high energies. 

Cosmic neutrinos may also cause air showers similar to cosmic rays, but deeper penetrating and thus distinguishable due to 
their large zenith angle [Q. At ~10^^ eV, such "horizontal" air showers are dominantly caused by atmospheric muons rather 
than cosmic neutrinos. Above 10^^ eV, however, the atmospheric background becomes low, and the air scintillation technique 
used, e.g., in the HiRes Fly's Eye detector or the Telescope Array ||74|], largely improves the detectability of horizontal air 
showers, providing much larger detection volumes than underground detectors. Above ^lO^^eV, the planned Pierre Auger 



Observatory 1 75 1 is expected to achieve considerable event rates, using the same technique in conjunction with ground arrays 
for particle detection [12||. The main caveat of the technique is the large atmospheric background — horizontal air showers 
produced by muon neutrinos can be easily confused with air showers from atmospheric muons generated by the prompt decay of 
charmed particles Horizontal air showers produced by electron neutrinos have the unique property to be mixed hadronic 

and electromagnetic showers [p4[|, which allows to determine distinctive triggering criteria for hybrid detectors like the Pierre 
Auger Observatory, reducing the background [p^. This makes electron neutrinos most interesting for UHE neutrino astronomy. 

The detection probability for deeply penetrating horizontal air showers can be expressed as the product of the hadronic 
neutrino cross section, a^p, and the detector acceptance for an horizontal air shower. Models predict correspondingly that 
CTjyp oc f for El, > 10^^ eV [[76[]. The detector acceptance for the Pierre Auger Observatory has been calculated as « 
for El, 10^^ eV, where horizontal air shower detection is expected to be more efficient than other techniques [p^]. This gives 
rise to /V dct E^^'^, which is the same dependence as in the case in water/ice Cherenkov experiments. For other air shower 
experiments, the dependence of the shower acceptance on energy might be different, but a slow rise in the UHE regime seems 
to be a common feature. 

The expressions for the energy dependence of the detection probability are highly approximate; exact results require expensive 
Monte-Carlo simulations, considering the detailed properties of the experiment. However, we note that the neutrino event rate per 
logarithmic energy interval, evaluated for the most common detector techniques, follows closely the neutrino power spectrum, 
4>v,dct/ ~ const. An exception is only the range between 10^^ cV and 10^^ cV, where horizontal air shower observations are 
still dominated by the atmospheric background and underground experiments affected by earth shadowing. 

APPENDIX D: TIME SCALES FOR FERMI ACCELERATION 

The time scale for first order Fermi acceleration at parallel shock fronts (defined as havin g th e flow direction parallel to the 
magnetic field lines, with perpendicular and obhque shocks defined accordingly) is given by pS] 
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iacc = -I -^(y- +xi3y+j , (Dl) 

where /3sh is the velocity of the shock in the comoving frame of the unshocked fluid, and and y+ are the ratios of diffusion 
coefficients parallel to the magnetic field, K^^, to the Bohm diffusion coefficient, viz., y — SK^^/r^^c < 1, taken in the regions 
upstream and downstream from the shock, respectively. The velocity jump at the shock in its comoving frame, xp — 
satisfies X/3 4 in nonrelativistic shocks (for an ideal gas with specific heat index of | ^^), and is x/3 = 3 in the ultra- 
relativistic limit [[77[]. Assuming y+ ^ = y, we find 6'f ~ PshV since Psh = P- by definition. For parallel shocks, there 
is no upper limit on the value of y, which can be interpreted as a measure for the strength of the turbulence in the magnetic 
field, y ~ (B/SB)^. In perpendicular shocks, one can show that 6'f ^ PshVli^ + J/^)' ^^'^ V limited from kinetic theory 



and isotropy requirements by j/ < [}~^ yTq]. In ultra-relativistic shocks, /3sh ~ 1, Eq. (Dl) is not strictly valid, but numerical 
simulations for both parallel and oblique relativistic shocks suggest 0^ ^ 1 |^], in correspondence with the result of Eq. ( ^T| ) 
for /3sh — y = ^- The time scale for second order Fermi acceleration is given by 



3cn, fSB 



(D2) 



where Va is the Alven speed in the plasma. Introducing y — {B/6B)^ as above, this yields tacc ^ yP^'^'t-L, with /3a = Wa/c. 
The condition for efficient scattering of the particles is /J^ <C y [^], which directly translates to 6'f ^ J/Z^a ^ ^ 1 for second 
order Fermi acceleration. 

Hence, tacc ~ = 27rrL/c, or 6'p ^ 1 gives a reasonable lower limit for the acceleration time scale of any kind of Fermi 
acceleration. It may be reached for acceleration in relativistic shock waves; in most cases, however, factors 0p > 10 would be 
more realistic. Throughout the paper, 6'f is treated as constant, i.e., which assumes that the diffusion coefficient is proportional 
to the Bohm diffusion coefficient. This is not true in other turbulence spectra, e.g., Kolmogorov turbulence, where 6'f 1 at 
the maximum energy im plie s (?f S> 1 at lower energies. While this can be important for the comparison of electron and proton 



acceleration time scales 1 80 1, it does not affect too much our results near the maximum proton energy. 



APPENDIX E: VHE/UHE NEUTRINO EMISSION FROM AGN CORES 

The first models which predicted considerable VHE/UHE neutrino fluxes from AGN assumed particle acceleration at shocks 
in the accretion flow onto the putative central black hole 1^,0,0. By assuming in these models that the power-law 2— lOkeV 
X-ray emission observed from AGN is produced by tt" -decay from p7 and pp interactions, and assuming that the observed 
diffuse X-ray background is entirely due to AGN, one can estimate the corresponding neutrino flux arising from pp and p^ 
interactions. While the assumption of shocks in the accretion flow near the black hole is more speculative than in jets (unlike the 
shocks in extended jets, the inner portions of AGN have never been imaged with sufficient angular resolution to infer shocks), 
core shocks are still the most cited class of models used to estimate expected event rates in neutrino experiments. Thus, it useful 
to investigate such models for self-consistency in the face of the pion and muon cooling effects discussed in this paper, which 
were not considered in the pubUshed results (see also more recent papers, e.g., [| 10|). 



In principle, we could incorporate these models into our general discussion, because the relevant sizes and time scales in 
AGN cores are also limited by variability. However, we will not write down here all the observational quantities, but use rather 
the physical parameters applied in the original papers. The only relevant quantity in the expression of the critical energy for 
pion and muon cooling (Eqs. (pO|)) is the value of the magnetic field, since AGN core models are obviously not Doppler boosted 
{V — 1). Assuming equipartition of the magnetic field and radiation energy densities, the magnetic field is generally taken 
as _B ^ [lO'^ G]£45^^^, where £45 is the UV luminosity of the AGN. (Note that the inverse dependence of the equipartition 
field on the luminosity arises from the fact that the luminosity scales with the distance of the shock from the Black Hole, 
and thus with the linear size of the acceleration region, as £ oc i?, leading to Uph oc £^^.) Applying Eqs.(|o|) this leads to 

£,^.syn{p) ^ [10^® eV]£45^, and £y syn(^) ^ [10^^ eV]£45^. Disregarding pion and muon cooling, the models of Stecker et al. 

predict a flat single-AGN neutrino power spectrum up to [2x 10^^ eV]£45^ followed by an exponential cutoff; Protheroe 
and Szabo [^|8|] find essentially the same result. The model of Sikora and Begelman predicts a sharp cutoff at about 10^^ eV 
due to their more conservative assumption for the acceleration time scale, 9^ ^ 100/3^^, rather than 9^ ~ f3~^ as assumed in 
the other models. The latter is the only model which is not modified by considering pion and muon cooling, while for the other 
models we see that the electron neutrino spectrum and 50% of the muon neutrino spectrum, arising from muon decay, steepen 
more than one order of magnitude lower in energy that previously assumed. 

For the prediction of detector event rates, the integrated diffuse neutrino background contributed by all AGN needs to be 
determined. Taking simple step functions as approximations for both the steepening induced by muon cooling and the exponen- 
tial cutoff induced by the maximum proton energy, and noting that the dependence of the maximum energy on the AGN X-ray 
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luminosity remains the same, we can derive a relation which allows to transform the result obtained disregarding muon cooling 
into the result expected when this effect is considered. The ratio of the integrated neutrino number flux, J-{£), to the unmodified 
single source spectrum, fo{S), for £ between the cutoff energies of the least and most luminous AGN, can thus be written as 



where PagnC-C) is the AGN luminosity function, and Oh (a;) is the Heaviside step function, 83(2;) = 1 for x > and Quix) — 
otherwise. Replacing £ ^ 2x 10^^ eV, as originally assumed, by £' — 10^^ eV as obtained from muon cooling, obviously leads 
to the relation Q'{£) = Q{££l£') w 0.05Q(£), since Q oc £-1 Since Q' / Q = T' IT, the diffuse flux of electron 

neutrinos (as well as muon neutrinos from muon decay) is reduced to about 5% of the value previously obtained, independent 
of the luminosity function used. Additionally, the exponential cutoff of the diffuse background, corresponding to the cutoff of 
the most luminous quasars, sets in already below 10^^ eV rather than at 10^^ eV. Similarly, the diffuse muon neutrino flux from 
pion decay is reduced to about 50% of the original value, so that the total VHE muon neutrino flux drops by about a factor of 3 
and cuts off at 10^* eV. 

Clearly, a more detailed calculation is required to obtain reliable flux rates under consideration of detailed spectral modification 
induced by pion and muon cooling, but our approximate results akeady show that the effect is important. In particular, we expect 
no considerable contribution from AGN cores to the electron neutrino spectrum in the energy range interesting for horizontal 
air shower measurements. No change of the predicted fluxes is expected in the energy range relevant for deep underwater or ice 
Cherenkov detectors, like Lake Baikal or AMANDA. In the interesting intermediate range, in particular relevant for the event 
prediction for proposed radio Cherenkov detectors, we obtain a moderately lower flux of muon neutrinos, and a severely reduced 
contribution of electron neutrinos. Since the model prediction are upper limits (constrained by the diffuse X-ray background), 
the drop in the rates cannot be balanced by adjusting astrophysical parameters. 
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